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ABSTRACT 

We present the preliminary results of a survey aimed at significantly increasing the range and 
completeness of the local census of spectroscopically confirmed white dwarfs. The current census 
of nearby white dwarfs is reasonably complete only to about 20 parsecs of the Sun, a volume that 
includes around 130 white dwarfs, a sample too small for detailed statistical analyses. This census is 
largely based on follow-up investigations of stars with very large proper motions. We describe here 
the basis of a method that will lead to a catalog of white dwarfs within 40 parsecs of the Sun and 
north of the celestial equator, thus increasing by a factor of 8 the extent of the northern sky census. 
White dwarf candidates are identified from the SUPERBLINK proper motion database, allowing us 
to investigate stars down to a proper motion limit /i > 40 mas yr -1 , while minimizing the kinematic 
bias for nearby objects. The selection criteria and distance estimates are based on a combination 
of color-magnitude and reduced proper motion diagrams. Our follow-up spectroscopic observation 
campaign has so far uncovered 193 new white dwarfs, among which we identify 127 DA (including 9 
DA+dM and 4 magnetic), 1 DB, 56 DC, 3 DQ, and 6 DZ stars. We perform a spectroscopic analysis 
on a subsample of 84 DAs, and provide their atmospheric parameters. In particular, we identify 11 
new white dwarfs with spectroscopic distances within 25 pc of the Sun, including 5 candidates to the 
D < 20 pc subset. 

Subject headings: Solar neighborhood - surveys - techniques: spectroscopic - white dwarfs - proper 
motions - stars: distances 



1. INTRODUCTION 

Statistics of the local white dwarf population, such as 
the space density, luminosity function, and mass distri- 
bution, are fundamental tools for understanding the evo- 
lution of the Galactic stellar populations and quantify - 



ing their ages (Oswalt et al. 1996 Leggett et al. 1998) 



Because of their low luminosities, obtaining a large and 
complete census of white dwarfs within a well-defined 
volume remains a challenge. The best volume that can 
be defined for a census of low-luminosity objects is the 
solar neighborhood, which alleviates the need for deep 
surveys, and also allows one to map out the sample in 
velocity space using readily available proper motions. 

A catalog and analysis of the sample o f white dwarfs 
within 20 pc of the Sun were presented by|Holberg et al.| 
and later refined by |Holberg et al.| ( |2008p and 
et all p009l . In light of these studies, the current 



(|2002l, 
Sion et 



census of nearby white dwarfs is believ ed to be 80% com - 



plete, and contains 127 white dwarfs (Sion et al. 2009 ). 
Every white dwarf suspected to lie within 20 parsecs of 
the Sun was analyzed in greater detail by |Giammichele| 



|et al. (2012), and 130 members ended up in their sample 
of local white dwarfs. Even if one assumes that the local 
sample is complete, the size of the sample is too small 
for detailed statistical analyses, and there is a need to ex- 
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tend the census and obtain a complete sample of white 
dwarfs fr om a larger volume. Suc h an effort was under- 



taken by Subasavage et al. ( 2009| and earlier references 



within) by measuring trigonometric parallaxes for new 
white dwarfs that are candidates of the 25 pc sample, 
as part of their DENSE project focused on o bjects in 



the southern hemisphere. Holberg et al. (2011) also an 
nounced that the complete sample of white dwarfs will 
be extended to 25 parsecs, thus doubling the volume of 
the local sample. Based on the space density of white 
dwarfs known within 10 pc of the Sun, [Subasavag e et al| 
( 2009 1 estimated that the census of white dwarfs within 
25 pc and with accurate trigonometric parallaxes is only 
~ 40% complete, and if we extend this horizon a little 
further — to 40 pc for instance — the census of white 
dwarfs remains largely incomplete. 

Nearby white dwarfs have been traditionally discovered 
in catalogs of stars with high proper motions. Major 
contribut i ons hav e been mad e, for instance , by |Luyten| 
( |1979a|b| , [Giclasj ( |1971[ ), and |Giclas et al| ( fl978| , who 



identified a significant number of faint, blue, high proper 
motion stars, and their pioneer work is still useful to to- 
day's astronomers. Indeed, in the first study dedicated 
to buildin g a complete census o f the local sample of white 



dwarfs by | Holberg et al 
form an important traction o 



LHS, G, and GD objects 
the 109 objects reported 
in that sample. Major contributions to the complete- 
ness of t he local white dwarf sam p le also come from th e 

( |2004l 



wor k of iVennes k Kawkal (12 0031), 

mo 



and |Kawka V cnncs (2006 ), w 

NLTT catalog oflSalim & Gouldl (|2003h , and in particu 



Kawka et al. 
surv eyed the revised 
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lar, identified eight new white dw arfs lying within 20 pc 
of the Sun. The contribution of Farihi et al. (2005) is 
also worth mentioning in this effort, as we ll as those of 



Subasavage et al.| ( |2007[ [2008 |2009), and |Sayres et al. 
(2012]), aimed at completing the 25 pc sample! 



But in order to extend the volume of our complete 
sample of white dwarfs, the first step is to identify 
nearby stars with smaller proper motions, the coolest 
(T e g ~ 3500 K) of which are extremely faint due to the 
intrinsic small radius of white dwarfs. With the goal 
to improve the statistics of the local white dwarf pop- 
ulation, we have been hunting for white dwarfs in the 
SUPERBLINK catalog. This catalog, which is based on 
a re-analysis of the Digitized Sky Surveys — with its 
20-45 yr baseline — is at least 95% complete for the en- 
tire northern sky down to V = 19.0, with a very low 
rate of spurious detection. It thus constitutes an ideal 
database from which to search for faint, high proper mo- 
tion objects such as nearby white dwarfs. Also, because 
of its low proper motion limit (fi > 0'.'04 yr -1 ), the SU- 
PERBLINK sample effectively eliminates the kinematic 
bias for stars in the immediate vicinity of the Sun, which 
is a known limitation of traditi onal catalogues such as 



the LHS catalog (/x > CK'5 yr 1 ; Luyten 



1979a} and the 

NLTT catalog (fj, > 0'.'18 yr" 1 ; |Luyten||1979bf . Hence, 
the SUPERBLINK catalog also represents a powerful 
tool for the study of the solar neighborhood. Search- 
ing this database should provide a complete sample of 
white dwarfs to a much larger distance limit. 

Also, the high completeness and deep magnitude limit 
of SUPERBLINK allows the detection of all white dwarfs 
down to the lu minosity function turn off, which occurs at 
L/Lq ~ 10~ 4 (Fontaine et al. 2001), up to a relatively 
large distance. For a 0.6 Mq white dwarf with a pure 
hydrogen atmosphere, for instance, this corresponds to 
T off = 5000 K, or M y =15.23. The limiting magnitude 
of V = 19 implies that SUPERBLINK should be de- 
tecting all white dwarfs down to the luminosity function 
turnoff to a distance of 56.7 pc from the Sun. The main 
question is what fraction of these stars are expected to 
have proper motions above the SUPERBLINK limit of 
[i > 0'.'04 yr -1 . Assuming that the distribution of ve- 
locities for white dwarfs to be the same as that of main- 
sequenc e stars in the vicin i ty of t he Sun, we can use Fig- 
ure 1 of Lepine k. Gaidos] pOll), which shows the kine- 
matic selection effects of SUPERBLINK by illustrating 
the fraction of stars in the Hipparcos catalog that would 
be selected with a proper motion cut of /i > 0'.'04 yr -1 
up to a given distance. At 56.7 pc, more than 90% of the 
stars are detected. This minimal kinematic bias there- 
fore allows one to detect most white dwarfs down to the 
luminosity function turnoff, and to perform a complete 
statistical analysis on a sample ~ 10 times larger than 
the current 20 pc census. 

The interest in the local population of white dwarf 
stars in not only statistical, but also astrophysical. In- 
deed, probing the solar neighborhood allows the detec- 
tion of faint, cool stars that would remain undetected at 
larger distances. Since the cool end of the white dwarf 
luminosity function is incomplete, obtaining a reliable 
estimate of the space density of white dwarf stars and 
comparing the luminosity function to models remains a 
challenge. The completion of the cool end of the white 



dwarf luminosity function would allow the accurate de- 
termination of the Galactic age and the verification of 
the white dwarf cooling t heory. Furthermore, many cool 
white dwarfs are peculiar ( Giammichele et al.|2012 ), and 
it is among them that we can expect to find transition 
objects that would allow us to establish the link between 
the different spectral types and to achieve a better un- 
derstandi ng of the white dwarf spectral evolution. The 
catalog of 
of stars o: 



Holberg et al. ( 2008 1 contains a large number 
particular astropnysical interest. For a de- 
tailed description of these stars, see | Giammichele et aL] 
( 2012[ ) and references therein. It is expected that surveys 
at 25 and 40 parsecs will unveil an even larger number 
of such key objects. 

In this paper, we search the SUPERBLINK catalog to 
extend significantly the census of wh i te dwa rfs in the so- 
lar neighborhood. Lepine & Shara (2005) have shown 
how reduced proper motion diagrams constructed from 
the SUPERBLINK catalog can produce a large num- 
ber of white dwarf candidates. We present here a more 
detailed search and identification of these white dwarfs 
through a large spectroscopic follow-up program. Our 
specific goal is to obtain spectral confirmation of all sus- 
pected white dwarfs within 40 parsecs of the Sun. Given 
the enormous amount of data and limited telescope ac- 
cess, we restrict ourselves to the northern part of the sky. 
The first step of this spectroscopic survey consists in the 
identification, observation, and classification of the white 
dwarf candidates. In Section 2, we present the catalog 
from which the candidates are obtained. We detail our 
selection method in Section 3, as well as distance esti- 
mates and candidate list in Section 4. Section 5 describes 
the results of our spectroscopic observation campaign, 
while a preliminary spectroscopic analysis, including the 
determination of atmospheric parameters, is provided in 
Section 6. Finally, a discussion follows in Section 7. A 
more thorough model atmosphere analysis of the atmo- 
spheric parameters of our complete survey of new white 
dwarfs within 40 parsecs will be reported in subsequent 
papers. 

2. PROPER MOTION AND PHOTOMETRIC DATABASE 

Our white dwarf candidates are identified from the SU- 
PERBLINK catalog of stars with proper motions fi > 40 
mas yr -1 . This catalog, based on a re-analysis of the 
Digitized Sky Surveys (which include POSS-I and POSS- 
II plate scans), is estimated to be > 95% complete in 
the northern hemisphere down to a visual magnitude 
of V = 19, but extends to V ~ 20 in many areas of 
higher Galactic latitudes. The current version of the 
catalog comprises 2,283,540 objects, all designated by 
the letters "PM I" followed by 10 characters based on 
the right ascension (a) and declination (6) of the ob- 
ject. The b asic search algorithm is described in |Lepine] 
et al. (2002), while quality control procedures, including 
cross-correlation with other catalogs and the compilation 
of astrom etric and photometric r esult s, are discussed at 
length in Lepine & Shara (|2005[ ) and Lepine & Gaidos 
(20111. A complete list of 61,977 northern sta rs with 
H > 150 mas yr -1 has already been published in Lepine 



fc Shara| ( |2005[ ) . We provide below a brief summary of 



the astrometric and photometric entries given in the cur- 
rent SUPERBLINK catalog. 
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Table 1 

Available Photometric Data 



Catalog 


Version 


Bands 


Counterparts 


Reference 


2MASS 
SDSS 


DR6 


JHK S 

ugriz 


1,472,665 
345,958 


Skrutskie et al.||2006) 


Adelman-McCarthy ct al. ( 2008 ) 


Hipparcos, Tycho-2 






118,000 


van ijccuwen (2UU7) 








[HWk ct al. (200U|| 


USNO-B1.0 




Bj, Rf,In 


1,567,461 


Monet ct al. (2U03) 


GALEX 


GR6 


FUV, NUV 


143,806 


Gil de Paz~et al.|l|2009| 



2.1. Astrometry 

SUPERBLINK provides coordinates on the Interna- 
tional Celestial Reference System for the 2000.0 epoch. 
For stars catalogued in Hipparcos, the positions are ex- 
trapolated to t he 20 00.0 epoch from the values given in 



van Leeuwen (2007), which are listed for the 1991.25 
epoch. Likewise, those not in Hipparcos but listed 
in Tycho-2 have their positions extrapolated from the 
proper motions listed i n Tycho-2, and if a star has a 



)03) 

"2M 



2003), its position 



counterpart in 2MASS ( Cutri et al 
is extrapolated from the position of the 2MASS counter- 
part. Finally, coordinates for stars without a Hipparcos 
or 2MASS counterparts are calculated by SUPERBLINK 
from the position of the stars on the POSS-II scans. The 
coordinates of those stars are thus less accurate but are 



gener ally within a few arcseconds (see Lepine & Shara 
2005 for details). 

SUPERBLINK also lists proper motions for each en- 
try, tabulated from three sources. When available 
proper motions are taken from the Hipparcos ca talog 
(|van L ccuwcn 2007} or from the Tycho-2 catalog (H0g 
et al. |2Tj00[ ). Otherwise, the proper motions listed are 



those measured in the SUPERBLINK proper motion sur- 
vey, based on the Digitized Sky Survey images. SU- 
PERBLINK ends up providing proper motions for more 
than 2 million stellar objects, and in particular, a to- 
tal of 1,567,461 stars with S > 0. From now on, when 
we mention the SUPERBLINK database, we refer to the 
northern part of the catalog. 

2.2. Photometric Data 

The construction of reduced proper motion diagrams 
requires, in addition to proper motion measurements, a 
set of photometric data in order to estimate the color 
of each star. Fortunately, the cross-correlation of SU- 
PERBLINK with other catalogs not only allows coor- 
dinates and proper motions to be measured with more 
accuracy, but it also provides a useful set of photometric 
data covering a large portion of the electromagnetic spec- 
trum. We describe these data in turn, and a summary is 
provided in Table [T] 

The Two Mi cron All Sky Survey (2MASS) Point 
Source Catalog ( |Skrutskie et al.|[2006 1 represents an ex- 
cellent source of near-infrared magnitudes for our targets 
in SUPERBLINK since the 2MASS survey covers the 
whole sky an d is complete down to J ~ 16.5. Lepine & 



[Shara] | 

white d 



2005| see their Figure 30) successfully showed that 
iwarts in SUPERBLINK could easily be separated 
from other stellar populations in a Hy vs V — J reduced 
proper motion diagram. For the present study, we used 
a version of the SUPERBLINK catalog in which 2MASS 
counterparts had already been found and assigned to 



1,472,665 of the stars (~ 94%), with the remainder hav- 
ing no detectable counterpart in 2MASS. These infrared 
J, H, and K$ magnitudes have a 0.02-0.03 mag accuracy 
down to 13th magnitude, and point sources are detected 
with S/N better than 10 fo r stars brighter than J = 15.9, 
H = 15.0, and K s = 14.3 dSkrutskie et al. 112006 1. 



The Sloan Digital Sky Survey (SDSS) also represents 
a useful source of photometric data, with ugriz photom- 
etry f rom the Data Release 6 ( Adelman-McCarthy et al. 
2008| ) for 345,958 counterparts in the SUPERBLINK cat- 
alog. The SDSS magnitudes have photometric uncer- 
tainties of roughly 1% in th e griz bands and 2% in u 



(Padmanabhan et al. 2008). This is by far the most 



accurate optical photometry available in our study, and 
will be especially useful to identify white dwarfs in the 
SUPERBLINK catalog. 

Optical photometry in the blue (Bt) and in the vi- 
sual (Vt) range are also extracted for 118,475 stars 
with counterparts in the Hipparcos and Tycho-2 cata- 
logs. Additional optical photo metry was also obt ained 
from the USNO-B1.0 database ( |Monet et aL||2003| ), pro- 
viding photographic magnitudes tor the totality of the 
catalog, i.e. 1,567,461 objects. However, for some en- 
tries, the photometry is available only for one or two 
bands. More specifically, Bj magnitudes are available 
for 1,390,471 objects, Rf magnitudes for 1,405,840, and 
In magnitudes for 912,550 objects. The blue Bj mag- 
nitudes are extracted mostly from scans of IllaJ plates 
from the Palomar Sky Surveys (POSS-I, POSS-II) and 
the Southern ESO Schmidt (SERC) Survey, the red R F 
magnitudes are extracted from scans of IIIaF plates 
from POSS-I and POSS-II and also from the Anglo- 
Australian Observatory red survey ( AAO-red) , while the 
near-infrared In magnitudes are extracted from the IVn 
plates from POSS-II and SERC. The B T and V T magni- 
tudes are more accurate (0.1 mag or better) than the pho- 
tographic magnitudes (typically 0.5 mag), but they are 
available only for the brightest stars in SUPERBLINK, 
while photographic magnitudes are available for every 
object. 

Finally, we also searched the sixth data release (G R6) 



of the GALEX database ( |Gil de Paz et al-1[2009|) and 
identified 147,096 counterparts to the SUPERBLINK 
objects (for a 5'.'0 search radius). The corresponding 
far-ultraviolet (FUV, 1350-1780 A) and near- ultraviolet 
(NUV, 1770-2730 A) magnitudes are particularly useful 
for the identification of blue objects, and in particular 
white dwarf stars. 



3. SELECTION OF THE CANDIDATES BASED ON 
REDUCED PROPER MOTION DIAGRAMS 
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Reduced proper motion diagrams (RPMD) are a par- 
ticularly efficient tool to identify white dwarf candidates 
with known proper motions (see , for instance, |Knox et al.| 
19991 lOppenheimer et aT1[200Tl |Vennes k KawkaPOO^ 
Carollo et al.|2006||Kilic et al.|2006l) . The reducedproper 
motion ot an object is denned as H m = m + 5 log /i + 5, 
where m is the apparent magnitude in some bandpass 
and fi is the proper motion measured in arcseconds per 
year. The reduced proper motion is analogous to the 
absolute magnitude M m = m + 51og7r + 5, where the 
trigonometric parallax 7r is replaced with the proper mo- 
tion fj, of the object. A reduced proper motion diagram 
is thus similar to a color-magnitude diagram, and white 
dwarfs occupy a similar location in the diagram, i.e. the 
bottom-left region. Furthermore, using the tangential ve- 
locity ftan = 4.74 /i7r -1 in units of km s _1 instead of the 
proper motion, we obtain H m = M m + 51ogw t an — 3.38, 
and each star population can be isolated based on the 
mean value of its tangential velocity. 

One major problem with the identification of white 
dwarf candidates using reduced proper motion diagrams 
is the contamination of the white dwarf region by other 
stellar pop ulations, and by h i gh-vel ocity subdwarfs in 
particular. Vennes & Kawka ( |2003[ ) showed, however, 
that this contamination can be substantially reduced by 
the inclusion of a criterion based on V — J. Similarly 
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Kilic et al. (2006) demonstrated that reduced proper mo 
tion diagrams are efficient for detecting cool white dwarfs 
only when the measured proper motions of all stellar 
populations are reliable, since subdwarfs with inaccu- 
rate proper motions can contaminate the other stellar 
populations, and notably the white dwarf region of the 
diagram. SUPERBLINK has an estimated false detec- 
tion level of less than 1% down to V — 19, but the false 
detection rate increases significantly for fainter sources. 
In our selection criteria, we thus restrict our search to 
that stars with V < 19. Fortunately, SUPERBLINK has 
a very high level of completeness for V < 19, exceed- 
ing 98% for most of the sky. We are thus confident that 
we can easily identify a significant fraction of the nearby 
white dwarfs using this technique. The next sections de- 
scribe the four reduced proper motion diagrams we used 
to identify white dwarf candidates in SUPERBLINK, in 
an effort to take advantage of the whole set of photomet- 
ric information available. The order in which these are 
presented follows the order of their estimated efficiency 
at isolating the white dwarf population, starting with the 
most efficient one. 

3.1. RPMD Using ugriz Photometry 

We present in Figure [I] the H g vs g — z reduced 
proper motion diagram constructed from the 345,958 
SUPERBLINK objects in the northern hemisphere with 
ugriz photometry available from the 6th Data Release 
of the Sloan Digital Sky Survey. As a result of the rela- 
tively high accuracy of the SDSS magnitudes, the white 
dwarf population is particularly well separated from the 
other populations in this diagram. 

To verify the accuracy of our procedure, we also display 
as red dots in Figure [I] the sample of white dwarfs taken 
from the 2008 May electronic version of the Catalogue 
of Spectroscopically Identified White Dwarfs 1 (McCook 
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Figure 1. Reduced proper motions diagram (H g ,g — z) for the 
345,958 stars in the SUPERBLINK catalog (p > 40 mas yr" 1 ) 
with counterparts in the SDSS 6th Data Release. The red dashed 
lines show the limits used to define the location of the white dwarfs. 
These limits include 4929 white dwarf candidates, i.e. 1.4% of the 
sample of SUPERBLINK stars with SDSS photometry. The spec- 
troscopically confirmed white dwarfs from the WD Catalog are 
shown as large red dots. 



1 http:/ /www. astronomy. villanova.cdu/WDCatalog/index. html 



k Sion| [19991 hereafter WD Catalog) with 5 > also 
found inSUPERBLINK and with g and z photometry 
available. We first note that 191 spectroscopically con- 
firmed white dwarfs lie in the expected region near the 
bottom left of the diagram, and that a small number of 
white dwarfs are color outliers. More precisely, 3 white 
dwarfs overlap the subdwarf region, and 4 are found 
in the redder, main sequence portion of the diagram. 
Two of th ese are actually b inary systems: 0855+604.1 
is a DBQ fcreenst ein 1969) and 0855+604.2 is a DCE? 
(Eggen k Grccnstcin 1965), while 1133+358 is an unre- 
solved DC+dM ( |Greenstein||l976 |. 

Fortunately, there are enough spectroscopically identi- 
fied white dwarfs that are well separated from the other 
populations to allow us to define selection criteria for the 
white dwarf area. These criteria are defined by the need 
to include as many white dwarf candidates as possible, 
while trying to keep the contamination from subdwarfs 
to a minimum. As a general criterion, the area occupied 
by the white dwarfs must include at least 80% of the 
spectroscopically confirmed white dwarfs. In the present 
case, this limit between the halo subdwarfs and the white 
dwarf area is defined by the following linear equation 
H g = 3.5(g — z) + 16.5. The slope and y-intercept are 
chosen in a very conservative manner, and include 96% 
of the WD Catalog sample with measured ugriz photom- 
etry. The reason why we recover almost all white dwarfs 
from the WD Catalog is that the known white dwarf pop- 
ulation is well separated in this reduced proper motion 
diagram, and only a few objects are color outliers. How- 
ever, in the upper part of the diagram (H g < 14) there is 
again some confusion between the different stellar popu- 
lations defined by our linear equation, so we simply ap- 
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ply an additional cutoff at B g = 14 based on the known 
white dwarf population, in order to keep the contamina- 
tion to a minimum. These adopted selection criteria for 
the H g vs g — z diagram are displayed in Figure [l] Out of 
the 345,958 SUPERBLINK stars with a counterpart in 
the 6th Data Release of the SDSS, about 4929 fall within 
the white dwarf selection limits. This number represents 
1.4% of the stars in the catalog with SDSS photometry, 
or 0.5% of the total number of objects in SUPERBLINK. 

3.2. RPMD Using GALEX Photometry 
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Figure 2. Reduced proper motions diagram (Hy, NUV — V) 
for the 147,096 stars in the SUPERBLINK catalog (/i > 40 mas 
yr _1 ) with counterparts in the 6th Data Release of the GALEX 
survey. The red dashed line shows the limits used to define the 
location of the white dwarfs. This limit includes 19,150 white 
dwarf candidates, i.e. 12.7% of the sample of SUPERBLINK stars 
with GALEX photometry. The spectroscopically confirmed white 
dwarfs from the WD Catalog are shown as large red dots. 



White dwarf stars are generally hotter than main se- 
quence or subdwarf stars, and since their atmospheres 
are usually devoid of heavy elements that could absorb 
the UV flux, they are moderately strong UV emitters and 
can easily be distinguished from non-degenerate stars in 
a reduced proper motion diagram built from GALEX 
photometry. We present in Figure [2] the By vs NUV — V 
diagram containing 147,096 stars in SUPERBLINK with 
NUV magnitudes measured by GALEX. NUV magni- 
tudes are used in this diagram since they are available 
for a much larger number of stars than the FUV magni- 
tudes. No corrections are applied for interstellar redden- 
ing, since according to the characterization of the Local 
Bubble of |Reis et al. (2011), the smallest distance to the 
wall of dust that causes extinction (E(b — y) > 0.040) is 
~ 80 pc. The reddening should thus not affect the white 
dwarf candidates of the local sample. 

Here and in the following diagrams, the V magnitudes 
are estimated from the relation 



as recommended by Lepine & Shara (2005), where Bj 
and Rf are photographic magnitudes taken from the 
USNO-B1.0 catalog. These estimated V magnitudes 
are believed to be accurate to ±0.5 mag ( |Lepine "fc 
Shara|20"05| . As mentioned earlier, while Vr magnitudes 
from the 'Lycho-2 catalog are more accurate than photo- 
graphic magnitudes, they are only available for a small 
number of SUPERBLINK objects, whereas Bj and Rp 
exist for the majority of our targets. Hence, despite the 
relatively large uncertainties of the V magnitudes em- 
ployed here, they have the advantage of being available 
for most of the entries in SUPERBLINK. 

As in Figure [T] the spectroscopically confirmed white 
dwarfs from the WD Catalog are also shown in red. Here, 
a total of 13 spectroscopically confirmed white dwarfs 
white dwarfs are scattered in areas normally occupied 
by other stellar populations. Most of these objects are 
cool degenerates, including 9 DA, 2 DC, 1 DZ (1705+030, 
Greenstein 1984) and 1 DQ star (1105+412, Koester and 
Knist, 2006). Most likely, these have very little UV flux 
and thus inaccurate NUV magnitudes. 

We define the slope and y-intercept of the line that 
characterizes the white dwarf region with the same cri- 
teria as before. Also, since the halo subdwarfs and main 
sequence stars are generally not as bright in the UV as 
white dwarfs are, there is no need to apply any further 
cutoff in By . Therefore, to be considered a white dwarf 
candidate, a star must have a reduced proper motion 
larger than By = (NUV - V) + 14. With this limit, 94% 
of the white dwarfs from the WD Catalog with measured 
NUV photometry are recovered. Finally, if this criterion 
is applied to the 147,096 SUPERBLINK objects with 
both NUV and V magnitudes, we obtain a sample of 
19,150 white dwarf candidates, which represent 12.7% of 
the stars with GALEX photometry in SUPERBLINK, or 
2.0% of the complete catalog. 

3.3. RPMD Using 2MASS Photometry 
As discussed above, the reduced proper motion dia- 



gram in By vs V — J, displayed in Figure 30 of |Lepine| 



& Shara (2005), showed that nearly 2000 white dwarfs 
could be identified in the SUPERBLINK catalog of stars 
with proper motions /i > 150 mas vr . Since the SU- 
PERBLINK catalog now includes stars with fi > 40 mas 
yr _1 , and because a significant fraction of its entries has a 
counterpart in the 2MASS catalog, such a diagram has an 
even greater potential for identifying white dwarf stars. 
The resulting By vs V — J diagram is shown in Figure 
[3j and contains 1,265,733 stars from SUPERBLINK with 
J magnitudes taken from 2MASS as well as Bj and Rp 
magnitudes from the USNO-B1.0 catalog. Here again, 
the V magnitudes are obtained from the empirical cali- 
bration given by Equation (1). 

The comparison with the spectroscopically confirmed 
white dwarfs from the WD Catalog reveals 20 objects 
that fall in a region to the right of that occupied by the 
bulk of white dwarfs. Among them, there are 9 DA, 4 
DC, 1 DQ, and 1 DZ star. This diagram also includes the 
largest number of multiple-star systems identified so far. 
Indeed, we f ind four WD +dM systems, a DB+DC binary 
(2058+342, piHnI|[2004|, and 0023+388, the triple star 
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Figure 3. Reduced proper motions diagram (Hy, V — J) for the 
1,265,733 stars in the SUPERBLINK catalog (fi > 40 mas yr" 1 ) 
with counterparts in the 2MASS catalog and for which an estima- 
tion of the V magnitude was possible. The red dashed lines show 
the limits used to define the location of the white dwarfs. These 
limits include 16,977 white dwarf candidates, i.e. 1.3% of the ini- 
tial sample. The spectroscopically confirmed white dwarfs from 
the WD Catalog are shown as large red dots. 
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Figure 4. Reduced proper motions diagram (Hy, V — In) for the 
878,847 stars in the SUPERBLINK catalog (fj. > 40 mas yr^ 1 ) with 
counterparts in the USNO-B1.0 catalog for all three Bj, Ftp and 
Im photographic magnitudes. The red dashed lines show the limits 
used to define the location of the white dwarfs. These limits include 
20,862 white dwarf candidates, i.e. 2.4% of the initial sample. The 
spectroscopically confirmed white dwarfs from the WD Catalog are 
shown as large red dots. 



system discussed earlier. 

Applying similar criteria as before, the white dwarf 
candidates are selected from Figure [3] if their reduced 
proper motion is larger than Hy = 3.5(U — J) + 13, 
with a cutoff of Hy > 14. These limits, displayed in 
Figure [3j recover 86% of the white dwarfs from the WD 
Catalog listed in SUPERBLINK and for which 2MASS 
and photographic magnitudes are available. This defines 
a sample of 16,977 white dwarf candidates out of the 
1,265,733 objects in the initial SUPERBLINK sample 
with 2MASS photometry, or a fraction of 1.3% (0.23% of 
the total initial sample). 

3.4. RPMD Using Photographic Magnitudes 

Even if USNO-B1.0 photographic magnitudes are less 
accurate than CCD photometry obtained in recent large 
surveys such as SDSS, they have the advantage of cov- 
ering the whole sky, and are thus available for all 1.6 
million SUPERBLINK objects. We present in Figure [4] 
the reduced proper motion diagram in Hy vs V — In 
constructed with V as defined by Equation (1). To be 
included in this diagram, photographic Bj, Rp, and In 
magnitudes must all b e available for each obje ct. How- 
ever, as discussed by Lepine & Shara (2005), not all 
USNO-B1.0 entries have magnitude information in all 
three bands. Whenever possible, Lepine & Shara tried 
to combine data if a USNO-B1.0 star appeared as more 
than one entry, but some sources remained without in- 
formation for one or more bands. As a result, Figure [4] 
includes 878, 847 stars out of the 1, 567, 461 entries in SU- 
PERBLINK with at least one photographic magnitude in 
USNO-B1.0. 

In Figure |1J the separation between the white dwarf 



and subdwarf populations is not as well defined as in 
the other reduced proper motion diagrams. The V and 
In filters are indeed too close in wavelength to allow an 
efficient separation of the two populations, as was the 
case with V — J or g — z, for instance. We also have to 
consider the fact that there is an uncertainty of ±0.5 mag 
in both V and In, which adds to the dispersion in U— In- 
The comparison of the white dwarfs in the WD Catalog 
displayed in Figure H] with those in Figure [3] using V — J 
reveals that most of the 55 outliers in V — In are also 
outliers in V— J, the difference between the two diagrams 
being that there is a larger concentration of outliers in 
V — In near the white dwarf locus. Once again, the 
measurement uncertainty in V— In is to blame. All of the 
interesting outliers have been discussed in the preceding 
sections and will not be repeated here. 

Despite this large contamination, we must still define 
some selection criteria, and as before, the slope and y- 
intercept of the limit between white dwarfs and halo sub- 
dwarfs is chosen to include 80% of the white dwarfs from 
the WD Catalog; a lower limit in Hy is also defined. Us- 
ing these criteria, the white dwarf candidates are selected 
from Figure [4] if their reduced proper motion is larger 
than H v = 3.5(V - I N ) + 13, with a cutoff of H v > 14. 
This time, since the contamination is much larger than 
in previous reduced proper motion diagrams, a more con- 
servative cut must be used, and no more than 80% of the 
white dwarfs from the WD Catalog are recovered. We 
finally end up with 20,862 white dwarf candidates out of 
the 878, 847 SUPERBLINK objects in the original sam- 
ple with USNO photometry, or a fraction of 2.2% of the 
total catalog. 

Combining the results from all four reduced proper mo- 
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tion diagrams, we obtain a total of 20,862 white dwarf 
candidates, since all of them have USNO photographic 
magnitudes, but only a fraction of them have data avail- 
able in other photometric systems. Some of these can- 
didates, however, can be found in up to four reduced 
proper motion diagrams. At this point, each candidate 
can be assigned an order of priority depending on the 
quality of the photometry used for its identification. 

3.5. A Priority Approach 

For each given star, data can be available for up to 
four photometric systems and their corresponding re- 
duced proper motion diagrams. Thus, it is possible that 
a star could be within the white dwarf region defined 
by our selection criteria in one diagram and outside in 
another diagram. We must therefore decide which photo- 
metric system should be prioritized to decide whether or 
not a star should be included in our final list of ~ 21, 000 
white dwarf candidates. For instance, ugriz magnitudes 
should take precedence over photographic magnitudes, 
since the former are much more accurate. In the follow- 
ing, we establish the order of priority for our four pho- 
tometric systems, based on their degree of photometric 
accuracy. 

SDSS magnitudes have a relatively high degree of ac- 
curacy and cover a wavelength range that allows an effi- 
cient separation of the white dwarfs from the other stel- 
lar populations, as discussed above. The reduced proper 
motion diagram based on ugriz magnitudes is arguably 
the most accurate, and it is therefore given the high- 
est priority. However, a comparison of our preliminary 
list of white dwarf candidates based on ugriz data with 
those found in the literatu re shows a significant contam- 
ination from subdwarfs. Kawka et al. (2004) success- 
fully reduce this contamination by including a color-color 
criterion to their reduced proper motion cut, namely 
(V — J) < 3.28( J - H) - 0.75. But such limits also get 
rid of the white dwarfs with red, main-sequence compan- 
ion s, which we would like to include. A l so, as observed 



by |Kilic et al.| ( |2006[ ) and |Sayres et ai] ( |2012[ ), any un- 
certainty in the proper motion measurements can lead to 
contamination. Also, spurious proper motions are pos- 
sible for stars near the faint magnitude limit of the SU- 
PERBLINK. With all these considerations in mind, we 
apply on to p of the criter i a in (H a ,g — z), a criterion 
inspired by |Kawka et al. (2004), but based on Figure 
[3j It is denned as follows: stars with V < 14 must 
also have Hy — 3.5(V — J) + 13. This limit on V re- 
moves the largest number of known main sequence stars 
from our list of candidates, while keeping all the known 
white dwarfs. Applying the criterion in (Hy, V — J) to 
fainter candidates results in the elimination of spectro- 
scopically confirmed white dwarfs, so it is not applied to 
those fainter stars. Finally, using the criteria defined in 
Section 3.1, combined with those defined in Section 3.3 
when V < 14, we retain a total of 4823 white dwarf can- 
didates. Unfortunately, the SDSS survey, at the epoch 
of the DR6, only covered about a quarter of the north- 
ern sky. Bright stars also tend to be saturated in SDSS, 
and we need to limit our selection to stars with u > 13 , 
g > 14, r > 14, i > 14, and z > 12 ( |York et al.||2000[ ). 
Consequently, there is no SDSS counterpart tor every star 
in the SUPERBLINK catalog, and we must therefore rely 
on other photometric systems to select our targets. 



In the absence of SDSS photometric data, GALEX UV 
photometry is used instead for the selection, whenever 
it is available. It is our second choice because of the 
corresponding photometric accuracy as well as the effi- 
ciency of the criteria in NUV— V to separate white dwarfs 
from other stellar populations. Here, we also apply the 
(Hy, V — J) criterion for stars with V < 14.0 in an effort 
to decontaminate our sample of candidates. The criteria 
and the method described in Section 3.2 led to the iden- 
tification of 8092 additional white dwarf candidates. The 
GALEX survey covers 80% of the sky, with special care 
taken to avoid the Galactic plane and Magellanic clouds, 
which could provide excess background flux in the UV 2 . 

In the absence of SDSS and GALEX data, 2MASS 
photometry combined with the criteria defined in Sec- 
tion 3.3 are used to identify 1132 additional white dwarf 
candidates with 6 > and V < 14.0. The 2MASS sur- 
vey is a precious source of photometric information since 
it covers practically the whole sky. However, since it is 
only complete down to J ~ 16.5, our fainter targets do 
not have a 2MASS counterpart. Finally, if there is no 
SDSS, GALEX, or 2MASS photometry available for a 
given target, we must rely on photographic magnitudes 
obtained from the USNO-B1.0 catalog. With the criteria 
defined in Section 3.4 and for V < 14.0, we identify an 
additional list of 6688 white dwarf candidates. 

All in all, a total of 20,735 white dwarf candidates are 
identified with the help of the four reduced proper mo- 
tion diagrams described in this section. This large num- 
ber of candidates amply justifies our decision to restrict 
our search to the northern hemisphere. Moreover, given 
our interest in establishing a census of white dwarfs in 
the solar neighborhood within 40 pc of the Sun, we must 
further restrain our list of candidates by evaluating pho- 
tometric distances for each object on our target list. 

4. A SAMPLE OF WHITE DWARF CANDIDATES WITHIN 
40 PC OF THE SUN 

4.1. Distances from Color- Magnitude Relations 

In the absence of trigonometric parallax measurements 
for most white dwarf candidates in our sample, we must 
rely on distances estimated using the only information 
available, which are the apparent magnitudes of each star 
in a set of bandpasses covering the ultraviolet to the near- 
infrared. In this section, we estimate photometric dis- 
tances from the distance modulus, m — M m — 5 log D — 5, 
where the absolute magnitude M m is determined from 
theoretical color-magnitude relations combined with a 
measured color index in some specified photometric sys- 
tem. To do so, we rely on synt hetic colors obtained 
using the procedure outlined in |Holberg fc Bergeron 
(|2006| based on the imp roved Vega fluxes taken from 
pohlin & Gilliland (2004). These color-magnitude rela- 



tions are available on our Web site 3 , or upon request on 
any photometric system. They cover a range of effec- 
tive temperature between T e g — 1500 K and 110,000 K 
and surface gravities between log g — 6.5 and 9.5 for 
hydrogen-rich atmospheres, and between T e g = 3500 K 
and 40,000 K and \ogg — 7.0 and 9.0 for helium-rich 
atmospheres. Note that we now include in our hydro- 
gen models the opacity from the red wing of Lya cal- 

2 http://galex.stsci.edu/GR6/ 

3 See http:/ /www. astro. umontreal.ca/~bergeron/CoolingModels 
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culated by Kowalski & Saumon (2006) and kindly pro- 
vided to us by P. Kowalski, which is known to affect 
the ffux in the ultraviolet region of the energy distribu- 
tion. The calculations of absolute magnitudes also re- 
quire the use of mass-radius relations for white dwarfs, 
which are b ased on evolutionary m odels similar to those 



described m 
g(He) = logM Ho 



Fontaine et al. (2001) but with C/O cores, 
/M* = 10 LV ~aiid g(H) 



10~ 2 and q(H) = 10~ 4 , which 
are representative of hydrogen-atmosphere white dwarfs, 
and <z(Hc) = 10 -2 and g(H) = 10~ 10 , which are rep- 
resentative of a helium-atmosphere. From these color- 
magnitude relations at constant log<? values, we can ob- 
tain the corresponding relations at constant mass values 
using the same evolutionary models as described above. 

We present below our results on the SDSS, GALEX, 
and 2MASS photometric systems, as well as on the 
USNO-B1.0 photographic system. Synthetic colors were 
calculated from the bandp asses available from th e SDSS 
website 4 and discussed in Fukugita et al 



( 1996t i while 



aL 
"fro: 



(2003) 



the 2MASS filters are described in Cohen et 
and the transmission functions were taken trom 

website 5 . Similarly, GALEX synthetic colors 



survey 

were calculated from the bandpas ses available from the 
GALEX websi t e 6 an d described in |Morrissey k GALEX 
Science Team (2004), and information about the Miters 



for the USNO-Bl.O are given in |Monet et aL] ( |2003| ) , and 
transmission curves are available Irom the Digitized Sky 
Survey website 7 . 
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4 http:/ /www.sdss.org/dr6/instruments/imager/#:nlters 

5 http:/ /www. ipac.caltech.edu/2mass/releases/allsky/doc/ 
sec6_4a.html 

6 http:/ /galexgi. gsfc.nasa.gov/docs/galcx/Documents/ 
PostLaunchResponseCurveData.html 

7 http:/ /www3.cadc-ccda.hia-iha.nrc-cnrc.gc.ca/dss/ 



Figure 5. Theoretical color-magnitude relations for pure 
hydrogen-atmosphere white dwarfs at 0.4, 0.6, and 0.8 Mq (dashed 
lines) and for pure helium-atmosphere white dwarfs at 0.6 Mq 
(dotted line). The pure hydrogen sequence at 0.6 Mq is used to 
determine absolute M g magnitudes for stars with SDSS photome- 
try. Also shown are the 50 white dwarfs from the WD Catalog with 
available ugriz photometry and trigonometric parallaxes from the 
Yale Parallax Catalog. 

We show in Figure [5] the theoretical M g vs g — z 
color-magnitude relation for hydrogen-atmosphere white 
dwarfs at 0.6 Mq together with similar sequences at 0.4 
Mq and 0.8 M Q , which are representati ve of the intrinsic 
mass distribution for DA stars (see, e.g., Gianninas et al. 
These are used below to evaluate the accuracy of 



[20TT | 



our color-magnitude calibration. Also shown in Figure 
[5] is a single 0.6 Mq helium- atmosphere sequence used 
to evaluate the influence of the unknown atmospheric 
composition on the color-magnitude relations. As can 
be seen, the 0.6 Mq helium sequence follows closely the 
corresponding hydrogen sequence in this particular di- 
agram, and it is thus perfectly justified to rely on the 
hydrogen-rich sequence only to evaluate the photometric 
distances to our objects. 

As an external verification of our color-magnitude re- 
lations, we also plot in Figure [5] the 50 spectroscopically 
confirmed white dwarfs from the WD Catalog that also 
have trigonometric paralla x measurements publish ed in 
the Yale Parallax Catalog ( |van Altena et aL 1995) and 
with SDSS photometry available; we consider here only 
stars with parallax uncertainties less than 30%. We also 
distinguish DA and non-DA stars. The absolute magni- 
tudes M g are directly obtained from the distance mod- 
ulus M g = g + 51og7r + 5. Our results show that the 
observed scatter with respect to the 0.6 Mq theoretical 
sequence is entirely consistent with that expected from 
the intrinsic white dwarf mass distribution, as indicated 
by the theoretical sequences at 0.4 and 0.8 Mq. This can 
be tested more quantitatively by measuring the mass of 
each star directly from the color-magnitude diagram us- 
ing our theoretical sequences and a simple interpolation 
scheme. We obtain a mean mass of (M) — 0.63 Mq, with 
a dispersion of a{M) — 0.20 Mq, entirely consi stent with 
the photome tric mass distribution obtained by Bergeron 



et al. (2001 see their section 5.3) for the sample of cool 
white dwarfs with trigonometric parallaxes from the Yale 
Parallax Catalog, (M) = 0.65 Mq, a(M) = 0.20 Mq. As 
also discussed by Bergeron et al., however, the dispersion 
expected from spectroscopic mass distributions are con- 
siderably smaller, typ ically a(M) ~ 0.15 Mq (see, e.g 



Gianninas et al.|2011 ) due to the increased sensitivity to 
log g ol the spectroscopic technique over the photometric 
method based on trigonometric parallax measurements. 
We thus conclude from this comparison that our M g vs 
g — z color-magnitude relation is well calibrated. 

Finally, we use the M g vs g — z theoretical relation for 
all white dwarf candidates in our sample with observed 
ugriz photometry to estimate a photometric distance for 
each object assuming a hydrogen-atmosphere and a mass 
of 0.6 Mq, the average mass for white dwarfs. The photo- 
metric distances for 4823 white dwarf candidates in our 
sample, identified from the (H g ,g — z) reduced proper 
motion diagram, are estimated in this manner. 

4.1.2. M v raNUV-F Calibration 
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Figure 6. Same as Figure [B] but for the determination of ab- 
solute My magnitudes for stars with GALEX photometry. Also 
shown are the 82 white dwarfs from the WD Catalog with avail- 
able GALEX photometry and trigonometric parallaxes from the 
Yale Parallax Catalog. The small dots correspond to USNO pho- 
tographic magnitudes while the larger symbols make use of Johnson 
V magnitudes. 



Mo, still consist ent with the values obtained by Berg 



(2001) 



eron et al 

We finally apply this calibration to all white dwarf 
candidates identified from the (£/y,NUV — V) reduced 
proper motion diagram, and estimate photometric dis- 
tances for 8092 stars with GALEX photometry (but with 
no SDSS counterparts). 
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In an approach similar to that described in the previ- 
ous section, we use the theoretical relation between the 
absolute magnitude My and the color index NUV — V to 
determine an absolute magnitude for every object with 
observed GALEX photometry, assuminga hydrogen at- 
mosphere at 0.6 Mq. We show in Figure [6] the theoretical 
My vs NUV — V color-magnitude relations for the same 
mass values and atmospheric compositions as above. In 
contrast with the results from the previous section, the 
helium sequence starts to differ from the hydrogen se- 
quence at 0.6 M Q for NUV - V > 2, or T cff < 8000 K, 
and becomes significantly different for NUV — V > 4, 
or T c g < 6000 K. However, since we do not expect to 
identify many white dwarfs cooler than ~ 6000 K on the 
basis of their UV magnitudes, it is justified to rely solely 
on the pure hydrogen sequence. 

Figure [6] also shows the 82 white dwarfs from the WD 
Catalog observed by GALEX, with trigonometric paral- 
lax measurements available in the Yale Parallax Catalog. 
Since we are mostly interested here in verifying the va- 
lidity of our color-magnitude relations, we want to use 
the best photometry available for each star in order to 
reduce the scatter related to the uncertainty of photo- 
graphic magnitudes. Hence, in addition to V magnitudes 
estimated from USNO photographic magnitudes (small 
dots in Figure [6]), we also show the results obtained using 
apparent V magnitudes measured on the Johnson pho- 
tometric system, and taken from the literature (larger 
symbols in Figure pi). This comparison is thus analogous 
to that shown in Figure [5j and we see once again that 
the bulk of white dwarfs is well contained between the 
0.4 Mq and 0.8 Mq theoretical sequences. Performing 
the same calculations as before, we find a mean mass 
of (M) = 0.67 Mq and a dispersion of <r(M) = 0.22 



V-J 

Figure 7. Same as Figure[5]but for the determination of absolute 
My magnitudes for stars with 2MASS photometry. Also shown are 
the 167 white dwarfs from the WD Catalog with available 2MASS 
photometry and trigonometric parallaxes from the Yale Parallax 
Catalog. The small dots correspond to USNO photographic magni- 
tudes while the larger symbols make use of Johnson V magnitudes. 



4.1.3. My vs V — J Calibration 

The theoretical relations between the absolute My 
magnitude and color index V — J are shown in Figure 
[71 Also displayed are the 167 white dwarfs from the WD 
(Jatalog for which a trigonometric parallax measurement 
in the Yale Parallax Catalog and a V — J color index 
are available (both Johnson V and USNO photographic 
magnitudes are displayed, as explained in the previous 
section). 

We note again in this figure that most of the points 
are contained between the 0.4 Mq and 0.8 Mq theoret- 
ical sequences, with (M) = 0.67 Mq and a{M) = 0.21 
Mq, and that the hydrogen- and helium-atmosphere se- 
quences agree sufficiently enough to assume a hydrogen- 
rich composition for all objects in our sample. We obtain 
the absolute magnitudes and photometric distances for 
all 1132 white dwarf candidates having counterparts in 
the 2MASS catalog (but with no SDSS or GALEX coun- 
terparts), assuming a hydrogen-atmosphere and a mass 
ofO.6M . 

4.1.4. My vs V — In Calibration 

Figure [8] displays our last color-magnitude relation, 
that between My and the photographic color index 
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Figure 8. Same as Figure[B]but for the determination of absolute 
My magnitudes for stars with USNO-B1.0 photographic magni- 
tudes. Also shown are the 151 white dwarfs from the WD Catalog 
with available USNO photographic magnitudes and trigonometric 
parallaxes from the Yale Parallax Catalog; the small dots corre- 
spond to USNO photographic magnitudes while the larger symbols 
make use of Johnson V magnitudes. The dot-dashed lines repre- 
sent the 0.5 mag dispersion around the pure hydrogen sequence 
at 0.6 Mq estimated from the accuracy of the Bj, Rp, and In 
magnitudes. 



V — In- Also shown are the 151 white dwarfs from the 
WD Catalog with trigonometric parallax measurements 
in the Yale Parallax Catalog and V — In color indices 
available. Again, both Johnson V and USNO photo- 
graphic magnitudes are used in this plot. Not unexpect- 
edly, the comparison between the observed My values 
and the theoretical sequences reveals a much larger scat- 
ter, larger than that expected from the intrinsic mass 
distribution alone. Indeed, we measure a mean mass of 
( M) = 0.58 Mp,, which i s 0.07 M Q lower than the value 



of Bergeron et al. (2001) for a similar sample, but more 
importantly, we find a significantly larger dispersion of 
cr(M) = 0.26 M Q , compared to the value of 0.20 M Q 
obtained by Bergeron et al. This is a direct consequence 
of the lesser accuracy of photographic magnitudes. In- 
deed, most of the scatter observed here is likely due to 
the 0.5 mag uncerta inty in photographic In magnitudes 
( Monet_et al. 2003). For instance, we also illustrate in 
Figure M the effect of a 0.5 mag error on the V — In 
color-index for the theoretical hydrogen sequence at 0.6 
Mq . Most of the points are then contained within these 
boundaries. The reliability of the color-magnitude rela- 
tion is therefore limited by the accuracy of the photo- 
graphic magnitudes. Unfortunately, these photographic 
magnitudes are the only information available to esti- 
mate photometric distances for the 6688 candidates with 
no SDSS, GALEX, or 2MASS counterparts. 

We finally conclude from the figures above that the 
color-magnitude relations derived from SDSS, GALEX, 
and 2MASS photometry are comparable in their level of 
accuracy, and that the least accurate photometric dis- 



tances, as one might expect, are those estimated from 
photographic magnitudes. 

4.2. Error on Photometric Distances 

There are 70 spectroscopically confirmed white dwarfs 
with measured trigonometric parallaxes recovered by our 
four reduced proper motion diagrams. For these 70 white 
dwarfs, we can estimate photometric distances using the 
photometric system corresponding to the reduced proper 
motion diagram where each object was identified. For 
instance, 19 objects were identified on the basis of their 
SDSS photometry, that is, they were selected from the 
(H g ,g — z) reduced proper motion diagram while their 
photometric distance was estimated using the M g vs 
(g — z) color-magnitude calibration. Similarly, 8 objects 
were identified with the help of GALEX photometry, 35 
from 2MASS photometry, and 8 from photographic mag- 
nitudes. We finally end up with a sample of 70 confirmed 
white dwarfs with distances measured from trigonomet- 
ric parallax measurements, where each star has also been 
identified in at least one of the reduced proper motion di- 
agrams, and with a corresponding photometric distance 
estimate. 
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Figure 9. Distances obtained from trigonometric parallaxes com- 
pared with photometric distances estimated from theoretical color- 
magnitude relations. The dots represent the 70 spectroscopically 
confirmed white dwarfs with parallax uncertainties less than 30% 
and also selected from our reduced proper motion diagrams. The 
solid line represents the 1:1 relation, while the dashed lines repre- 
sent the lcr dispersion of 8.5 pc resulting from the combined errors 
of the 4 color-magnitude relations. The dotted line indicates the 
+15 pc error adopted in our analysis to ensure all white dwarf can- 
didates within 40 parsecs of the Sun are included in our sample. 



The distances obtained from parallax measurements 
are compared to photometric distances in Figure [9j The 
dashed lines represent the average 8.5 pc (rms) disper- 
sion relative to the 1:1 relation estimated using the white 
dwarfs displayed in Figures[5]to[8]for all four sets of color- 
magnitude relations. Given the observed dispersion in 
Figure [9] and the fact that this dispersion appears to in- 
crease with distance, we choose to include in our list of 
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white dwarf candidates within 40 parsecs of the Sun all 
objects with an estimated photometric distance less than 
55 pc. We believe that this conservative buffer of 15 pc 
(dotted line in Figure [9| is enough to include all white 
dwarfs that could potentially lie within 40 pc of the Sun. 
Of course, we keep in mind that the main purpose of 
this calibration is purely to identify the nearest objects. 
Subsequent spectroscopic analyses are expected to pro- 
vide more accurate distances, and lead to an independent 
estimation of the error on the photometric distances. 

4.3. List of White Dwarf Candidates 

We are now able to determine photometric distances 
for each object on our list of white dwarf candidates fol- 
lowing the method described in the previous sections. If 
the estimated distance places it within 55 parsecs of the 
Sun, the object becomes part of our list of candidates for 
follow-up spectroscopy. This process leaves us with a list 
of 1978 spectroscopic targets. The sample can be fur- 
ther reduced by eliminating all previously known white 
dwarfs. To do so, the coordinates of each candidate are 
compared to those listed in the WD Catalog, and then 
with every star in the Simbad Astronomical Database 8 
within a 1 arcmin radius. This way, we found 499 white 
dwarfs in our candidate list that were previously known, 
and an additional 35 objects with a known spectral type 
that identifies them either as a main sequence stars or 
a background galaxy. The presence of these 9 galaxies 
in our candidate list indicates that there are apparently 
some spurious objects with false proper motions in the 
SUPERBLLNK catalog. From all these objects with a 
known spectral type, we can evaluate the contamination 
of our sample of white dwarf candidates to be ~ 8%. 

With the objects with a known spectral type removed 
from our candidate list, we finally end up with a list 
of 1341 targets for follow-up spectroscopy. This sam- 
ple divides into 268 candidates identified on the basis 
of SDSS photometry, 130 from GALEX photometry, 731 
from 2MASS photometry, and 212 from USNO-B1.0 pho- 
tographic magnitudes. The candidates identified with 
SDSS and GALEX colors are given first priority for 
follow-up spectroscopy, followed by stars with 2MASS 
photometry. Finally, objects with only photographic 
magnitudes available are given the lowest priority. 

The 268 white dwarf candidates with available SDSS 
photometry are shown in Figure 10 in a {u—g, g—r) color- 
color diagram, together with theoretical predictions for 
our pure hydrogen, pure helium, and DQ models. This 
two-color diagram reveals that our sample is dominated 
by cool white dwarfs. We also expect that an impor- 
tant fraction of these cool objects will have a hydrogen- 
dominated atmosphere, and that some candidates are 
most likely DQ stars. 

The 1341 white dwarf candidates identified in SU- 
PERBLINK are also displayed in Figure [II] The upper 
panel shows their space distribution using a cylindrical 
equal-area projection of the equatorial coordinates, while 
in the lower panel, the same distribution is binned into a 
histogram to illustrate more clearly the space density as 
a function of right ascension. This figure shows that our 
white dwarf candidates are distributed uniformly on the 
sky, without any drop in surface density near the galac- 

8 http://simweb.u-strasbg.fr/ 
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Figure 10. (u — g, g — r) color-color diagram showing the 268 
white dwarf candidates with ugriz photometry available. The 76 
white dwarfs spectroscopically confirmed in Section 5 are shown 
with various symbols explained in the legend, while those without 
spectroscopic data are shown with red dots. The solid curves rep- 
resent pure hydrogen model atmospheres at logg = 7.0, 8.0, and 
9.0 (from bottom to top); effective temperatures are indicated in 
units of 10 3 K. The dashed curve corresponds to pure helium atmo- 
spheres at logg = 8.0, and the dotted lines represent DQ models 
for 5 different compositions, from logC/He = —9.0 to —5.0. 

tic plane. The spectroscopically confirmed white dwarfs 
from the WD Catalog are also shown for comparison. 
Given this characteristic of our survey, we hope to in- 
crease significantly the completeness of the local white 
dwarf sample in this region, which has often been ne- 
glected in the past due to field crowding near the Galactic 
plane. 

5. SPECTROSCOPIC FOLLOW-UP OF THE CANDIDATES 

We provide in this section a status report of our on- 
going spectroscopic survey, and present results for a first 
set of 422 objects from our target list, or nearly a third 
of our complete sample. 

5.1. Spectroscopic Observations 

Optical spectra have been obtained with the Steward 
Observatory 2.3-m telescope, and the NOAO Mayall 4-m 
and 2.1-m telescopes, during 8 different observing runs 
between 2009 May and 2010 October. The adopted con- 
figurations allow a spectral coverage of AA3200-5300 and 
AA3800-6700, at an intermediate resolution of ~ 6 A 
FWHM. Spectra were first obtained at low signal-to- 
noise ratio (S/N ~ 25), which is sufficient to identify 
main sequence objects, but also represents the lower limit 
required to obtain reliable model fits to the spectral lines. 
Table[2]summarizes our spectroscopic observing runs and 
observational setups. 

As a result of our spectroscopic observations, 193 
newly identified white dwarfs have been spectroscopi- 
cally confirmed. Among these, 68 were identified on the 
basis of SDSS photometry, 18 from GALEX, 70 from 
2MASS, and 12 from USNO photographic magnitudes. 
The remaining 25 objects were discovered using an 
earlier version of our selection method, based on criteria 
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Figure 11. Upper panel: Equal cylindrical projection of the equa- 
torial coordinates for the sample of 1341 white dwarf candidates 
identified from SUPERBLINK (solid circles) compared with the 
sample of 499 stars from the WD Catalog recovered by our selec- 
tion criteria (open circles). Also shown by the bold solid line is 
the region of the galactic plane. Lower panel: Space density as a 
function of right ascension, normalized to the total number of stars 
in each bin (both lines are thus on a comparable scale). The solid 
line represents the 1341 white dwarf candidates, while the dashed 
line corresponds to the white dwarfs from the WD Catalog. 

different from those adopted in the final version. These 
objects have revised photometric distance estimates 
beyond 55 pc, and are thus not included in our final list 
of 1341 candidates. 



5.2. Spectroscopic Content 

Our spectroscopic follow-up observations identify 193 
new white dwarfs, listed in Tables [3 and [4] Table [3] pro- 
vides astrometric data as well as NLTT and SDSs des- 
ignations, when available, while Table [2] lists the same 



objects but with the available photometry and adopted 
spectral types. In the earlier presentation of our results 
(Limoges et al. 20101, the ' LSPM J' notation was used 
tor the object names, as in Lepine fc Shara| ( [2005| , but 
the the proper motion limit for stars in SUPERBLINK 
has lowered to /x > 0.04"yr _1 , and the entries were ho- 
mogenized with a 'PM F designation. In summary, this 
subsample contains 127 DA (among which 9 DA+dM 
and 4 magnetic), 1 DB, 56 DC, 3 DQ, and 6 DZ white 
dwarfs. 

Once we have a confirmed white dwarf, it is possible to 
improve upon our preliminary distance estimates, which 
were based on approximate V magnitudes and color- 
magnitude relations, by making use of the full set of pho- 
tometric data. Here we rely o n the so-called photometric 



method described at length in Giammichele et al. (2012 



and references therein) where the available photometry 
for each star is fitted with theoretical fluxes, properly 
averaged over each bandpass. Both T cff and the solid 
angle ir(R/D) 2 are considered free parameters, where R 
is the radius of the star and D its distance from Earth. 
We assume a value of log 5 = 8 and corresponding radius 
R, and obtain directly the distance D. These improved 
photometric distances are given in Table [4] 

When this survey was undertaken in 2009 May, none of 
these objects had a white dwarf classification or a spec- 
tral type available in the literature. But since then, white 
dwarf i dentifications have be e n reported in |Kilic et al. 



2010 , Vennes et al. 



ayres et al -12012) 



Tonry et al. 



IKilic et al.] 
( |2012p , and 



( 2011D , ^ 

We identified these stars in Table 

tand chose to leave these objects in our sample since 
ey have been discovered independently. We also want 
to point out that even though Table [3] contains 25 stars 
with existing SDSS spectra, all our targets have been 
identified using our own spectroscopic data. 

Optical spectra from our subsample of DA stars cov- 
ering H8 t o H/3 — or Ha when available — are shown in 



Figure 12 Note that 14106+0245 (right panel of Figure 



12 d, second object from the top) is a DAZ, and that our 



subsample of 4 magnetic DA white dwarfs (05158+2839, 
06019+3726, 06513+6242, and 15164+2803) are dis- 



played in Figure 
white dwarfs wit 



13] Our survey also detected 9 new DA 
a an M dwarf companion; these are plot- 
ted separately in Figure [14) This was quite unexpected, 
since the cuts in the reduced proper motion diagrams 
were chosen in order to avoid main sequence stars. As 
a consequence, we avoided all objects that are bright in 
the infrared portion of the spectrum. As discussed ear- 
lier, however, this is true for our criteria in g ~ z, V — J, 
and V — In, but not for our criteria in NUV — V, which 
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Wavelength (A) 



Figure 12. (a) - Optical spectra for our sample of DA white dwarfs from SUPERBLINK. The spectra are displayed in order of decreasing 
equivalent width of H/3, from upper left to bottom right, and shifted vertically for clarity. The Ha line is also shown, when available, and 
normalized to a continuum set to unity. 



is efficient for detecting blue objects. And indeed, all 9 
DA+dM systems were detected using this last reduced 
proper motion diagram. 

The DB, DQ, and DZ stars in our sample are plot- 
ted together in Figure [15] Unfortunately, the observa- 
tional setup with the NOAO telescopes does not allow 
the coverage of wavelengths shorter than ~ 3900 A, while 
covering Ha simultaneously. Calcium lines can still be 
easily identified, however, but additional spectroscopic 
observations near the ~ 3700 A region are currently be- 
ing obtained in order to perform a proper model atmo- 
sphere analysis of these DZ stars. In the case of DQ 
white dwarfs, one particular object, 12476+0646, ex- 



hibits the pre ssure-shifted car bon lines characteristic of 
DQpec stars ( |Kowalski 2010[ ). Finally, we display our 
featureless DC spectra m Figure 16 in order of their Right 
Ascension. 

Some of the results of our spectroscopic observations 
are summarized in the color-color diagram shown in Fig- 
ure |10[ where we identify the various spectral types of 
the 76 confirmed white dwarfs with available SDSS col- 
ors. As expected, the DQ stars are located in the appro- 
priate region of the (u — g, g — r) diagram, and in the 
next phase of the survey, we plan to use this characteris- 
tic to identify all possible DQ stars in SUPERBLINK. 
We also note the presence of a DA + dM system at 
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Figure 12. (b) - continued. 

(u — g,g — r) ~ (2.8, 1.2), in the redder part of the di- 
agram. Finally, the sample of DC stars follows the the- 
oretical, pure hydrogen sequence, with only one outlier 
near u — g — 2.0, giving us a preliminary indication of 
the atmospheric composition even before performing a 
full analysis of their energy distribution. 

In the following section, we present a preliminary spec- 
troscopic analysis of the DA component of our survey. 

6. ATMOSPHERIC PARAMETER DETERMINATION OF DA 

STARS 

The coolest white dwarfs in our sample are either fea- 
tureless, or present too few spectral lines for a proper 
spectroscopic analysis, and the determination of their at- 



mospheric parameters (T e ff, logg) can only be achieved 
from an analy sis of their photometr ic energy distribu- 
tion (see, e.g., Bergeron et al. 1997). At the moment, 
not enough photometric information is available to pro- 
ceed with a homogeneous analysis of the coolest objects 
in our sample, and we are still securing the appropriate 
optical and infrared photometry for cool DA, DC, DQ, 
and DZ stars, the results of which will be reported in 
subsequent papers. We thus restrict our determination 
of the atmospheric parameters to the subsample of 84 
spectroscopically confirmed DA stars for which the spec- 
troscopic technique can be successfully applied. 
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Figure 13. Our subsamplc of magnetic DA white dwarfs, shifted 
vertically for clarity. 



X 
3 



o 
> 



o 
K 



05492+5747 




23229+3358 




6.1. Theoretical Framework 

Our model atmospheres and synthetic spectra for DA 
stars are built fro m the model atm o sphere code orig- 
inally described in jBergeron et al. (1995) and refer- 
ences therein, with re c ent im provements discussed in 
Tremblay & Bergeron (2009). These are pure hydro- 
gen, plane-parallel model atmospheres, with non-local 
thermodynamic equilibrium effects explicitly taken into 
account above T e s — 30, 000 K, and energy transport 
by convection is included in cooler models following the 
ML2/a = 0.7 prescription of the mixing-length theory. 
The theoretical s pectra are calculated within the occupa- 
tion formalism of Hummer & Mihalas ( 1988 ), which pro- 
vides a detailed treatment of the level populations as well 
as a consistent description of bound-bound and bound- 
free opacities. We also rely on the improved c alculations 
for the Stark broad ening of hydrogen lines from |Tremblay| 
fc Bergeron|( |2009| ), which include nonideal perturbations 
from protons and electrons directly inside the line profile 
calculations. Our model grid covers a range of effective 
temperature between T e s = 1500 K and 120,000 K, and 
log g values between 6.0 and 9.5. 

Our fitt ing technique is based on the approach pio 
neere d by IBergeron et al.l (119921 see also ILiebert et al. 



2005 1 , which relies on the nonlinear least-squares method 



of Levenberg-Marquardt ( Press et al.|1986 |. The optical 
spectrum of each star, as well as the model spectra (con- 
volved with a Gaussian instrumental profile), are first 
normalized to a continuum set to unity. The calculation 
of x 2 is then carried out in terms of these normalized line 
profiles only. The atmospheric parameters - T c g, log g - 
are considered free parameters in the fitting procedure. 

Special care needs to be taken in the case of DA stars 
with an unresolved M dwarf companion in order to re- 



Wavelength (A) 

Figure 14. Our sample of binary systems composed of a DA white 
dwarf and a M dwarf companion. The spectra are shifted vertically 
for clarity. The Ha and H(B line cores of 04032+2520E (not to be 
confused with the DA 04032+2520W), 06556+5920, 16171+0530, 
and 23283+3319 are contaminated by emission from the M dwarf. 

duce the contamination of the white dwarf spectrum by 
the companion. When the contamination affects only 
H/3, and sometimes H7 as well, we simply exclude these 
lines from the fit (e.g., 05280+4820 and 23229+3358). 
At other times, emission lines from the M dwarf are also 
observed in the center of the Balmer lines, in which case 
the line centers are also simply excluded from our fitting 
procedure (e.g., 04586+6209). A similar approach was 
also adopted if the flux contribution from the M dwarf 
is too important and "fills up" the Balmer line cores, re- 
sulting in predicted lines that are too deep (06556+5920, 
16171+0530, and 23283+3319). In some cases, however, 
the white dwarf spectrum is too contaminated by the 
M dwarf companion to be fitted with the simple ap- 
proach described above (e.g., 04032+2325E — not to 
be confused with the DA star 04032+2325W — and 
11036+1555 9 ), and a more robust fitti ng procedure us- 
ing M dwarf templates will be required ( Gianninas et al. 
2011 1. These results will be presented elsewhere. 



6.2. Spectroscopic Results 

Even though the spectroscopic technique is arguably 
the most accurate method for measuring the atmospheric 
parameters of DA stars, it has an important drawback 

9 In the case of 11036+1555, we even detect the 4226 A line from 
the M dwarf in the white dwarf spectrum. 



1G 



81551+4103 [ 



> 



0) 




02476+4138 



22418+0432 



Wavelength (A) 



Figure 15. Spectra of white dwarfs in our sample with helium- 
rich atmospheres. The panel on the left displays the DZ stars, 
while the panel on the right shows our only DB star, as well as 3 
DQ stars. 



at low effective temperatures (T e g < 13, 000 K) where 
spectroscopic values of log g are significantly larger than 
thos e of hotter DA stars, t he so-called high- log g p roblem 
(see Trem blay et al. | 2010| and references therein). Trem 
bfay et ai.| ( |2011bp showed that this high- log g problem 
is actually related to the limitations of the mixing-length 
theory used to describe the convective energy transport 
in DA stars, and that more realistic, 3D hydrodynami- 
cal model atmospheres are required in order to obtain a 
surface gravity distribution that resembles that of hot- 
ter radiative-atmosphere DA stars. Since these spurious 
high- log q values aff e ct dir ectly the estimated distances, 
Giammichele et al. (2012) derived an empirical procc- 
dure (see their Section 5 and Figure 16) to correct the 
log<? values based on the DA stars in the Data Release 



4 of the Sloan Dig ital Sky Survey, analyzed by |Trem- 
We adopt a similar approach her 



blay et al. 


(2011a) 


and apply t 


leir log 



T eS = 7000 K and 14,000 K. 

The spectroscopic fits for our subsample of 84 DA stars 
are displayed in Figure |17| The corresponding atmo- 

fheric parameters (T c g and log g) are reported in Table 
together with the stellar mass (M/M ), absolute ab- 
solute visual magnitude (My), luminosity (L/L@), esti- 
mated visual magnitude (V), spectroscopic distance (D), 
and white dwarf cooling time (r). Whenever necessary, 
we rely on the same evolutionary models as those de- 
scribed above to derive these quantities. In principle, 
the spectroscopic distance can be obtained directly from 
the distance modulus, by combining the theoretical ab- 
solute magnitude in a single given bandpass with the 
observed magnitude in the same bandpass. However, 



since the photometric errors can be large in some sys- 
tems we used — the USNO photographic magnitudes in 
particular — , we estimated the spectroscopic distances 
by using the full set of photometry available for each star, 
and calculated an average spectroscopic distance, prop- 
erly weighted by the photometric uncertainties in each 
bandpass. This is equivalent to using the photometric 
method described above but by forcing the effective tem- 
perature at the spectroscopic value, thus fitting only the 
solid angle ir(R/D) 2 , where R is the radius of the star 
determined from the spectroscopic logg value. In doing 
so, we also fold in the uncertainty of the spectroscopic 
logg measurement. 

The mass distribution for the DA stars in our sample 
is displayed in Figure [18] as a function of effective tem- 
perature. This figure clearly illustrates the efficiency of 
our survey to identify white dwarfs using reduced proper 
motion diagrams even at very low effective temperatures. 
We also distinguish with various color codes the criteria 
used in our survey to discover each white dwarf, allow- 
ing us to study the impact of one particular photometric 
system on the selection process as a function of temper- 
ature. For instance, white dwarfs identified on the basis 
of GALEX photometry extend down to relatively low ef- 
fective temperatures. Indeed, the observed photometric 
sequence allows us to apply our selection criteria down 
to NUV - V = 6.5 (see Figure [Hf, or T eS - 5300 K. 
Similarly, white dwarfs identified on the basis of ugriz 
photometry are mostly found at the low end of the tem- 
perature distribution. Most SDSS targets are intrinsi- 
cally faint, and thus include an impressive amount of cool 
white dwarfs that can only only be identified through the 
use of reduced proper motion diagrams. Surprisingly, 
white dwarfs identified on the basis of 2MASS photome- 
try are found at all temperatures. This is due to the fact 
that our photometric sequences allow us to apply our 
color criteria as blue as V — J = —0.5 (see Figure u\, or 
T c ff ~ 20, 000 K. Finally, only a few white dwarfs in this 
subsample were identified on the basis of USNO photo- 
graphic magnitudes. From these results, we can conclude 
that even though SDSS represents the most reliable pho- 
tometric data set, GALEX, 2MASS, and even photomet- 
ric magnitudes are also required to identify white dwarfs 
over the complete range of effective temperature. 

The mass distribution of DA white dwarfs in our sub- 
sample, regardless of their temperature, is displayed in 
Figure [T9] The mean mass of these 84 DA stars is 
0.689 M Q with a standard deviation of a = 0.172 M fl , a 
value significantly large r than the value obtained by |Gi- 



ammichele et al. ( 2012 1 for the DA white dwarfs witinn 
20 pc of the Sun (HM7 M with a = 0.171 M ). One 
obvious difference is that we do not include here the white 
dwarfs already known in the literature. Most likely these 
are brighter, intrinsically more luminous, and probably 
have larger radii and thus lower masses. The mass distri- 
bution of the 37 DA white dwarfs within 40 pc of the Sun 
displayed in Figure 19 (shaded histogram) actually shows 
an important high- mass component (see also Figure 18). 
These high-mass white dwarfs, with their small stellar 
radii, are intrinsically less luminous than their normal- 
mass counterparts, and they are thus more abundant in 
a volume- limited sample, such as the local neighborhood, 
than in a magnitude-limited sample. Our results indicate 
that we are successfully recovering these high-mass white 
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Figure 16. Our sample of featureless DC stars. All spectra are normalized to a continuum set to unity and are offset from each other by 
a factor of 0.9. 



dwarfs in our surv ey, often missing in magnitude- limited 



2005 in the case of the 



surveys (see, e.g., |Liebert et aL 
PG survey) 



7. DISCUSSION 

7.1. Comparison of Spectroscopic and Photometric 
Distances 

During the target selection process, distances were es- 
timated using approximate V magnitudes together with 
various color-magnitude relations, displayed in Figures [5] 
to [8] These distance estimates were later improved by 
comparing theoretical average fluxes to the set of avail- 
able photometry, properly weighted by their uncertainty. 
At that point, we simply assumed a surface gravity of 
log g — 8.0, and considered both T c r and the solid an- 
gle as free parameters. These estimates are referred to 
as photometric distances. The spectroscopic analysis, on 
the other hand, provides spectroscopic distances, where 
for a given star, theoretical absolute magnitudes are com- 
puted from the spectroscopic values of T c g and log g, and 
compared to the set of available photometry. In both 
cases, if only the V photographic magnitude is available, 
the 0.5 magnitude error will introduce a 23% uncertainty 
on the estimated spectroscopic distance. If additional 



photometry is available, however, this distance uncer- 
tainty can be significantly reduced (see Table [5]). 

The comparison between photometric ana spectro- 
scopic distances for the DA white dwarfs in our sample 
is displayed in Figure 20 We restrict this comparison to 
T c ff > 7000 K since the B aimer lines in cooler objects 
become too weak to be fitted properly with the spectro- 
scopic method, yielding spurious log g values at low tem- 
peratures (see Figure 18) and corresponding distances. 
The dotted lines in Figure [20| represent a ±23% difference 
between both estimates (i.e., the maximum error on spec- 
troscopic distances obtained from photographic magni- 
tudes, as discussed in the previous paragraph). The bulk 
of stars is generally found within these limits. Part of the 
observed dispersion in Figure [20| can be attributed to the 
intrinsic mass distribution of our sample. Indeed, all our 
color-magnitude calibrations assumed a typical mass of 
0.6 Mq, but as shown in Figures [5] to |8| there is an in- 
trinsic dispersion in absolute magnitude due to the mass 
(or radius) distribution of white dwarfs. In particular, 
white dwarfs with very high (M > 0.9 M ) spectro- 
scopic masses yield photometric distances that are over- 
estimated; these are identified with a different symbol in 
Figure [20} 
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Figure 17. (a) - Fits to the optical spectra of the DA stars in our sample. The lines range from H/3 (bottom) to H8 (top), each offset 
vertically by a factor of 0.2. Theoretical line profiles shown in green are not used in the fitting procedure. 



Another source of scatter is due to the presence of M 
dwarf companions, which make the system brighter at vi- 
sual and infrared magnitudes compared with single DA 
stars. Since these magnitudes were used to estimate their 
photometric distance (My vs V — J), this can easily ac- 
count for the large discrepancies with spectroscopic dis- 
tances. Indeed, in the spectroscopic distance calculation, 
the less accurate magnitudes weigh less, and the more 
accurate JHK photometry dominates the distance solu- 
tion. Finally, as noted in Table [5j we have certain doubts 
about the cross-correlation with the GALEX database 
for a handful of stars in our sample. For these objects, 
the GALEX photometry is inconsistent with the rest of 
the spectral energy distribution, and they had to be omit- 
ted from the fits used to estimate the spectroscopic dis- 
tances. However, as in the previous cases, these colors 
were used to obtain our initial distance es tima te. 

To summarize, most objects in Figure [20] are found 



between the ±23% dispersion in distance, and the stars 
falling outside these limits can be separated into three 
categories: DA stars with M dwarf companions, high- 
mass white dwarfs, and stars with large photometric un- 
certainties (see the corresponding error bars in Figure 



20). The previously estimated 15 pc error is thus enough 
to identify white dwarfs with reasonably accurate pho- 
tometry, and we thus conclude that searching at 55 pc 
in order to find all white dwarfs within 40 pc is realistic, 
especially when photometry such as SDSS, GALEX, or 
2MASS is available. 

Our preliminary spectroscopic analysis of DA stars 
presented in Table [5] yields 11 white dwarf candidates 
within 25 pc of the sun, including 5 candidates within 
the 20 pc sample. Incidentally, a few of these objects 
already have a parallax measurement avail able. Indeed, 



21134 +0727 (G25-20) has a parallax from |Dahn et al. 
(|1988|, 7T = 0'.'0411±0 / .'0038 yr" 1 , placing it within 25 pc. 
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Figure 17. (b) -continued. 

Also, if 22118+5649 is a common proper motion compan- 
ion to LTT 16500, as it is suspected to be (Subasavage et 
al. 2012, private communication), then it has a parallax 
of 7T = 0"02677 ± 0'.'00018 yr" 1 (or D = 37.4 pc), and 
thus not a member of the 20 pc sample, while still within 
40 pc of the Sun. Finally, a private communication from 
J. Subasavage confirms that 16325+0851 is indeed within 
25 pc of the Sun. So even though spectroscopic distances 
are more accurate than the previous photometric esti- 
mates, the only way to confirm the membership of white 
dwarfs to the local sample is through trigonometric par- 
allax measurements. Such measurements would not only 
provide reliable distances, but would also yield mass de- 
terminations for the coolest objects in our sample ana- 
lyzed with the photometric technique. 

7.2. Success Rate of Discovery 

The absolute visual magnitudes (estimated from the 
calculated V magnitudes and photometric distances) for 



the 193 new white dwarfs identified i n o ur survey are 
plotted in the upper panel of Figure 21 as a function 
of photometric distance. Also shown are the 499 white 
dwarfs in SUPERBLINK already known in the literature. 
The candidates still without spectroscopic confirmation 
are displayed separately in the lower panel; the objects 
selected on the basis of their USNO photographic mag- 
nitudes are considered second priority targets because of 
their higher probability of being contaminants from the 
main sequence. In each panel, the dashed lines repre- 
sent lines of constant apparent V magnitude. We note 
that the white dwarfs identified in our survey are dom- 
inated by objects fainter than V — 16, and that most 
of them are found at photometric distances larger than 
20 pc. This is not surprising since the census of white 
dwarfs within 20 pc of the Sun is believ ed to be at least 
90% complete flGiammichele et al.||2012[ ). There are still 
a few white dwarf candidates on our target list within 
20 pc that have no spectroscopic data, due to observa- 
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Figure 17. (c) - continued 

tional constraints, but these stars are currently on our 
high priority list. 

From the results shown in the upper panel of Figure 
|21| we can determine that the ratio of new to known 
white dwarfs is 193/499 ~ 39%. Also, out of the 286 
candidates observed, 220 are confirmed white dwarfs (27 
in the literature 10 and 193 in our survey), for a success 
rate of 77%. This number is close to the 80% expected 
from our selection criteria, and we conclude that our sur- 
vey is quite efficient for recovering the missing fraction 
of white dwarfs in the solar n eigh borhood. 

The lower panel of Figure [2l] reveals that a signifi- 
cant fraction of our remaining white dwarf candidates are 
fainter than V = 17 (590 objects fainter versus 329 ob- 
jects brighter). The spectroscopic identification of these 

10 Some spectra of spectroscopically confirmed white dwarfs were 
secured by us and will be used in our next paper as part of a study 
of the total white dwarf content of SUPERBLINK within 40 pc of 
the Sun. 



stars with 2 to 4-m telescopes requires integration times 
on the order of an hour under excellent weather con- 
ditions. The candidates deserving spectroscopic follow- 
up must then be carefully chosen, and our high-priority 
list now includes 89 of these faint candidates and 186 
"bright" targets, for a total of 275 high-priority targets, 
excluding 120 objects we already observed after 2010 Oc- 
tober and that are still being reduced. Future observa- 
tions will be dedicated to the follow-up of these high- 
priority white dwarf candidates, in particular those iden- 
tified on the basis of SDSS or GALEX photometry. 

7.3. Increasing the completeness of the current census 

We have already established the success rate of our 
spectroscopic survey, and we are now interested in its 
completeness. First of all, our white dwarf sample is di- 
rectly affected by the completeness of the SUPERBLINK 
catalog, which is high because of its low proper motion 
limit (// > 0'.'04 yr _1 ) which minimizes the kinematics 
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Figure 17. (c) - continued 

bias. To illustrate the effect of proper motion on kine- 
matics, we plot in Figure [22] the transverse motions Vt 
(i.e. the projected motions on the plane of the sky, where 
Vt = 4.47/id) for all stars in our sample as a function of 



the photometric distance D, as ca l culate d in Section 4.1. 
As explained in Lepine & Gaidos (2011), a star at 50 pc 
from the Sun with /z > CK'04 yr _i has a transverse veloc- 
ity v t < 9.48 km s" 1 , which will occur with a probability 
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Figure 18. Mass as a function of effective temperature for a sub- 
sample of 84 DA white dwarfs with spectroscopic mass determina- 
tions. All stars are identified with a different color based on the 
photometric system from which they were discovered. The dotted 
lines represent the theoretical isochrones for our C/O core evolu- 
tionary models with thick hydrogen layers, corresponding to the 
white dwarf cooling age in units of Gyr. The dashed line indicates 
the median mass of DA white dwarfs, 0.594 Mq, as determined by 
|Tremblay et al.| | |2011a) . 
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Figure 19. Mass distribution for the DA white dwarfs in our 
sample. The solid line histogram shows the distribution for the 
84 DA stars with spectroscopic masses available, while the shaded 
histogram corresponds to the subsample of 37 objects with spec- 
troscopic distances less than 40 pc from the Sun. The mean values 
and standard deviations are given in the figure. 

of about 10% for stars in the solar neighborhood (see 
their Section 2.2 and Figure 1). Their diagram shows 
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Figure 20. Comparison of photometric and spectroscopic dis- 
tances for white dwarfs with T e g > 7000 K, as defined in the text. 
The solid line represents the 1:1 correspondence, while the dotted 
lines represent a ±23% dispersion. 
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Figure 21. Absolute magnitude as a function of photometric dis- 
tance. In the upper panel, the filled circles represent the 193 new 
white dwarfs identified in our survey, while the open circles cor- 
respond to the 499 white dwarfs already known in the literature. 
Also shown by dashed lines in the figure are lines of constant ap- 
parent V magnitudes. The white dwarf candidates in our survey 
without spectroscopic confirmation are shown in the lower panel. 
The lower-priority candidates (those identified on the basis USNO 
photographic magnitudes) are shown with cross symbols. 

that operating with a proper motion limit of /i > Cf!l50 
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Figure 22. Transverse velocity as a function of distance, show- 
ing the kinematic bias due to the proper motion limit of the SU- 
PERB LINK catalog (fi > 40 mas yr — 1 ). In the upper panel, the 
filled circles represent the 193 new white dwarfs identified in our 
survey, while the open circles correspond to the 499 white dwarfs 
already known in the literature. Also shown by dashed lines in 
the figure are lines of constant apparent proper motion (in units 
of arcsec yr — 1 ). The white dwarf candidates in our survey with- 
out spectroscopic confirmation are shown in the lower panel. The 
lower-priority candidates (those identified on the basis USNO pho- 
tographic magnitudes) are shown with cross symbols. 



yr" 1 will only detect half of the stars at 40 pc and very 
few stars (only those with very large components of mo- 
tion) at 100 pc. However, a sample with a proper motion 
limit fx > 0'.'04 yr -1 will include ~ 95% of the stars at 
40 pc and ~ 70% of the stars at 100 pc. 

Hence, in terms of new white dwarf identification as 
a function of proper motion, we find that f or u > 0'.'5 
yr" 1 (see corresponding dashed line in Figure 22 1, which 
corresponds to the limit of the LHS survey, the ratio 
of new to known white dwarfs is 8.8%, while this ratio 
reaches 43.5% for 0'.'5 > [i > 0'.'18 yr" 1 (where the lower 
proper motion limit is that of the NLTT survey), and it 
then drops slightly to 41.3% for 0'.'18 < fi < 0'.'04 yr" 1 . 
Our survey is thus more efficient for proper motions lower 
than the LHS limit, but our results also demonstrate that 
the sample of white dwarfs with /i > 0'.'5 yr" 1 could host 
up to 7 more white dwarfs. Previous searches for white 
dwarfs within the NLTT limit were also incomplete, since 
33 of our new identifications have an NLTT designation. 
Note also that the NLTT and LHS appear to be > 80% 
complete down to the 19th magnitude, but only in the 
Completed Palomar Region (CPR), i.e. for S > -32.5° 



and outside a band ±10° of the Galactic plane (Lepine 
fc Shara|[2005| 



Spectroscopic distances were obtained for 84 out of the 
193 newly identified white dwarfs, while a preliminary 
photometric analysis (not presented here) of a subsample 
of the coolest objects was performed for another 78 white 
dwarfs. From this combined analysis of 162 white dwarfs, 
we find that 126 objects are within 55 pc of the Sun, and 
93 within 40 pc. The spectroscopic analysis of 1151 DA 



stars by Gianninas et al. ( 2011 ) contains 223 white dwarfs 
from the WD Catalog located in the northern hemisphere 
whose spectroscopic distances are within 55 pc from the 
Sun, and 121 within 40 pc. Using this latter survey, our 
ratio of new to known white dwarfs is estimated at 56% 
within 55 pc, and 77% within 40 pc. This difference 
in ratio will most likely be reduced when all candidates 
between 40 and 55 pc are observed. It was also mentioned 
in Section 7.2 that the ratio of new to known white dwarfs 
from the literature was 39%, while our success rate in 
detecting white dwarfs (both new and known) is 77%. 
Our survey is thus efficient for recovering white dwarfs 
that are already found in the literature as well as new 
identifications. 

In spite of the success of our survey, the first sample 
of newly identified white dwarfs presented in this paper 
is far from complete, but the survey has not reached its 
limit yet. Our analysis represents the first results of an 
ongoing effort, and more data are still being collected 
and analyzed. Moreover, SUPERBLINK is currently be- 
ing cross-correlated with the SDSS DR7 and GALEX 
GR7, providing additional high-quality photometric in- 
formation to replace USNO magnitudes in our selection 
process. This will eventually result in more high-priority 
candidates, and will also help in the identification of DQ 
and DZ stars, which separate well in color-color dia- 
grams, as we showed earlier, but only when such color 
information is available. Finally, these new magnitudes 
may also complete the set of photometry for white dwarfs 
in the SUPERBLINK catalog, allowing fits to the energy 
distribution of the cool white dwarfs that cannot be an- 
alyzed spectroscopically. Our future catalog will provide 
more candidates for parallax measurements, as well as 
more cool, massive, magnetic, and astrophysically chal- 
lenging white dwarfs, while being at least 80% complete 
within 40 pc of the Sun. We will also be able to provide 
statistics on the Solar Neighborhood based on a sample 
of white dwarfs large enough to reduce the uncertainties 
related to small number statistics. 
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14277+0532 
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+28:03:22. 


.2 





,107 


-0.106 


-0.010 


2 


15206+3903 




J152038.34+39034 


15 


:20:38 


.36 


+39:03:49. 


.3 





,169 


-0.095 


0.139 




15263+2936 




J152621.08+29362 


15 


:26:21 


.06 


+29:36:22. 







,188 


0.086 


0.167 





Table 3 — Continued 









RA 


DEC 




fiRA 


/iDE 




PM I 


NLTT 


OIJOO 


(<^200o) 


(^200o) 


(II -1\ 

v y r ) 


in -i\ 

v y r ) 


///,, r -i\ 




15342+0218 




J153417.49+02184 


15:34:17.50 


+02:18:48.0 


0.143 


-0.118 


-0.080 


2,4 


15359+2125 




J153554.02+21250 


15:35:54.06 


+21:25:06.6 


0.248 


-0.142 


0.203 




15494+4802 




J154927. 19+48022 


15:49:27.20 


+48:02:29.0 


0.176 


0.108 


-0.139 


2,4 


15589+0417 




J155857.68+04170 


15:58:57.69 


+04:17:05.1 


0.071 


-0.025 


-0.066 


2 


16053+5556 




J160520.60+55564 


16:05:20.69 


+55:56:45.4 


0.274 


-0.222 


0.160 




16096+4735 




J160939.78+47350 


16:09:39.74 


+47:35:10.7 


0.062 


0.043 


-0.044 




16171+0530 




J161710.56+05303 


16:17:10.56 


+05:30:38.2 


0.097 


0.011 


-0.096 




16264+1938 




J162626. 28+19384 


16:26:26.33 


+19:38:39.3 


0.275 


—0.174 


0.213 




16325+0851 




J163233.27+08512 


16:32:33.17 


+08:51:22.6 


0.375 


0.271 


-0.259 


4 


16335+5231 


43101 


J163332. 24+52314 


16:33:32.14 


+52:31:50.0 


0.265 


0.103 


-0.244 




16477+2636 




J164744.70+26364 


16:47:44.74 


+26:36:46.3 


0.204 


-0.151 


0.137 


2 


16546+5742 






16:54:41.01 


+57:42:16.3 


0.066 


-0.058 


0.031 




17027+1022 




J170246. 29+10224 


17:02:46.29 


+10:22:41.2 


0.255 


0.077 


0.243 




17052+0423 






17:05:13.25 


+04:23:45.3 


0.168 


0.009 


-0.167 




17238+0458 






17:23:49.57 


+04:58:49.1 


0.149 


-0.041 


-0.143 




17283+0211 






17:28:18.47 


+02:11:10.2 


0.224 


-0.029 


-0.223 




17335+7949 






17:33:34.29 


+79:49:16.3 


0.397 


0.394 


-0.045 




17417+2401 






17:41:46.23 


+24:01:49.6 


0.360 


— 0.030 


—0.359 




17433+1434S 






17:43:22.97 


+14:34:52.4 


0.386 


-0.252 


0.293 




17471+2859 






17:47:08.30 


+28:59:09.8 


0.198 


-0.121 


-0.157 




18014+5049 




J180127.50+50495 


18:01:27.49 


+50:49:58.2 


0.288 


-0.108 


0.267 




18073+0357 






18:07:23.37 


+03:57:01.7 


0.121 


-0.012 


-0.121 




18138+2119 






18:13:48.55 


+21:19:20.6 


0.050 


0.019 


0.046 




18199+1739 






18:19:59.23 


+17:39:19.0 


0.230 


0.216 


0.080 




18572+2026 






18:57:13.95 


+20:26:28.4 


0.316 


-0.013 


-0.316 




19128+5343 






19:12:48.56 


+53:43:13.5 


0.205 


0.148 


0.142 




19132+2949 






19:13:16.53 


+29:49:28.2 


0.287 


-0.056 


0.281 




19146+1428 






19:14:36.11 


+14:28:25.4 


0.422 


—0.114 


—0.406 




19464+0937 






19:46:29.80 


+09:37:14.9 


0.574 


0.060 


-0.571 




19493+0747 






19:49:23.51 


+07:47:31.7 


0.184 


0.137 


-0.122 




20062+5902 




J200616.93+59022 


20:06:16.96 


+59:02:27.3 


0.108 


-0.019 


-0.106 




20069+6143 




J200655.03+61430 


20:06:54.89 


+61:43:10.3 


0.238 


0.168 


-0.168 




20223+8333 






20:22:21.50 


+83:33:55.6 


0.208 


—0.127 


0.166 




20235+7001 






20:23:29.87 


+70:01:55.5 


0.320 


-0.260 


-0.186 




20597+5517 




J205945. 11+55173 


20:59:44.91 


+55:17:29.9 


0.498 


0.274 


0.417 




21077+0740 






21:07:45.05 


+07:40:44.3 


0.116 


-0.019 


-0.114 




21117+0120 






21:11:46.37 


+01:20:54.4 


0.068 


0.058 


-0.035 




21134+0727 






21:13:28.92 


+07:27:04.2 


0.341 


0.321 


0.115 




21222+0413 




J212212.31+04135 


21:22:12.35 


+04:13:56.7 


0.425 


-0.089 


-0.415 


4 


21384+1856 




J213827.28+18564 


21:38:27.28 


+18:56:41.9 


0.169 


-0.141 


-0.093 




21420+2252 






21:42:03.90 


+22:52:29.0 


0.186 


0.118 


-0.144 




21429+0805 






21:42:54.49 


+08:05:27.4 


0.218 


-0.098 


-0.195 




21492+0415 






21:49:13.60 


+04:15:50.2 


0.227 


0.011 


—0.226 




21551+4103 






21:55:06.36 


+41:03:06.6 


0.193 


0.172 


0.087 




21597+2936 






21:59:46.94 


+29:36:41.7 


0.097 


0.097 


0.005 


1 


22105+4532 






22:10:34.53 


+45:32:40.0 


0.361 


-0.122 


-0.340 




22118+5649 






22:11:53.64 


+56:49:46.4 


0.250 


0.249 


0.022 




22276+1753 


53908 




22:27:40.36 


+17:53:21.4 


0.219 


0.203 


—0.081 




22299+3024 






22:29:58.03 


+30:24:10.6 


0.122 


0.050 


-0.111 




22331+2610 






22:33:11.59 


+26:10:14.9 


0.089 


0.085 


0.026 




22418+0432 






22:41:53.18 


+04:33:00.2 


0.226 


0.115 


-0.194 




22447+1513W 






22:44:43.02 


+15:13:46.3 


0.291 


0.162 


0.241 




22497+3623 


54984 




22:49:46.43 


+36:23:22.8 


0.370 


0.330 


—0.167 




22595+5717 






22:59:33.81 


+57:17:57.3 


0.334 


-0.316 


-0.109 




23003+2204 






23:00:21.96 


+22:04:16.0 


0.225 


0.201 


0.102 




23027+4312 






23:02:44.89 


+43:12:47.7 


0.181 


-0.174 


-0.051 




23056+4334 






23:05:39.06 


+43:34:03.8 


0.270 


-0.120 


-0.230 




23098+5506E 






23:09:58.53 


+55:06:49.1 


0.408 


0.406 


0.040 




23160+0559 




J231605.60+05594 


23:16:05.46 


+05:59:46.6 


0.178 


0.175 


0.034 




23229+3358 






23:22:59.82 


+33:58:47.1 


0.179 


0.128 


0.125 




23234+7255 






23:23:28.24 


+72:55:07.4 


0.088 


0.069 


0.055 




23243+2835 






23:24:18.85 


+28:35:55.5 


0.134 


-0.050 


-0.124 




23283+3319 






23:28:18.03 


+33:19:31.1 


0.147 


0.053 


-0.137 




23389+2101E 






23:38:56.29 


+21:01:18.2 


0.338 


0.290 


0.174 
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RA 


DEC 


Mtot 




iuDE 


PM I 


NLTT SDSS 


(J2000) 


(J2000) 


(V- 1 ) 


CV- 1 ) 


("yr -1 ) Notes 


23390+5316 




23:39:00.48 


+53:16:00.0 


0.467 


0.466 


0.033 


23462+1158 




23:46:12.50 


+11:58:49.7 


0.136 


-0.107 


-0.084 


23475+0304 


J234735. 20+03043 


23:47:35.11 


+03:04:31.9 


0.278 


0.177 


0.215 


23478+0223 


J234753.64+02234 


23:47:53.76 


+02:23:40.8 


0.224 


-0.183 


-0.129 



Note. — (1) Also in Vermes et al. (2011); (2) Classified using our own data, but SDSS spectrum exists; (3) Also in Kilic ct al. (2010); (4) Also 
in Sayrcs ct al. (2012); (5) Also in Tonry ct al. (2012). 
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Table 5 

Atmospheric Parameters of DA White Dwarfs from SUPERBLINK 



PM I T cB (K) log g M/Mq M v log L/L G V a D(pc) log r Notes 



00331+4742S 
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(0.05) 
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(0.03) 


11 


14 
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96 


14 


75 


55 
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6 


(1 


1) 


9.15 


00532+3927 


7520 


( 44) 


8.00 


(0.09) 


0.60 


(0.05) 


13 


26 


-3 


33 


17 


05 


45 


2 


(1 


3) 


9.10 


01043+3816 


16380 


(102) 


8.05 


(0.05) 


0.64 


(0.03) 


11 


19 


-1 


99 


15 


35 


58 


3 


(1 


6) 


8.21 


01382+4442 


17180 


(157) 


7.95 


(0.05) 


0.59 


(0.03) 


10 


96 


-1 


85 


15 


57 


69 


3 


(2 


0) 


8.06 


01486+3615 


6470 


(138) 


8.65 


(0.28) 


1.01 


(0.17) 


14 


87 


-4 


01 


16 


93 


24 


2 


(° 


9) 


9.68 


01565+2955 


15170 


(121) 


8.00 


(0.05) 


0.62 


(0.03) 


11 


25 


-2 


10 


14 


64 


48 


3 


(1 


3) 


8.28 


02230+5544 


6930 


( 49) 


8.38 


(0.10) 


0.83 


(0.07) 


14 


13 


-3 


70 


16 


76 


35 


8 


(1 


1) 


9.53 


02557+2106S 


6410 


( 93) 


8.44 


(0.20) 


0.87 


(0.13) 


14 


54 


-3 


88 


17 


25 


37 


5 


(1 


2) 


9.64 


02562+4954 


6370 


(166) 


9.02 


(0.31) 


1.20 


(0.13) 


15 


65 


-4 


33 


17 


14 


16 


7 


(° 


7) 


9.65 


03127+2218 


6620 


( 63) 


7.61 


(0.16) 


0.38 


(0.07) 


13 


25 


-3 


34 


16 


99 


56 


9 


(1 


5) 


9.03 


03433+1958 


7160 


( 47) 


8.62 


(0.09) 


0.99 


(0.06) 


14 


41 


-3 


81 


15 


80 


18 


4 


(° 


7) 


9.60 


03473+4358 


13580 


(290) 


7.91 


(0.05) 


0.56 


(0.03) 


11 


32 


-2 


24 


14 


16 


33 


9 


(1 


0) 


8.37 


04032+2520W 


22590 


(246) 


8.11 


(0.05) 


0.69 


(0.03) 


10 


72 


-1 


46 


14 


89 


98 


1 


(2 


9) 


7.72 


04339+2827 


13860 


(318) 


8.01 


(0.05) 


0.61 


(0.03) 


11 


42 


-2 


26 


15 


75 


75 


9 


(2 


2) 


8.40 


04389+6351 


13050 


(626) 


7.67 


(0.10) 


0.44 


(0.05) 


11 


06 


-2 


17 


16 


85 


83 


5 


(2 


3) 


8.28 


04523+2519 


22180 


(195) 


8.55 


(0.05) 


0.96 


(0.03) 


11 


48 


-1 


79 


14 


87 


47 


7 


(1 


6) 


8.20 


04586+6209 


11610 


(967) 


8.83 


(0.12) 


1.12 


(0.06) 


13 


14 


-3 


13 


Id 


71 


24 


9 


(° 


9) 


9.23 


05025+5401 


11410 


( 70) 


8.01 


(0.05) 


0.61 


(0.03) 


11 


79 


-2 


60 


15 


19 


46 


9 


(1 


4) 


8.64 


05280+4850 


11120 


(380) 


8.44 


(0.14) 


0.88 


(0.09) 


12 


53 


-2 


92 


16 


04 


75 





(2 


4) 


8.95 


05492+5747 


13510 


(410) 


8.04 


(0.07) 


0.63 


(0.04) 


11 


52 


-2 


33 


15 


15 


75 


5 


(2 


1) 


8.46 


06026+0904 


6020 


( 42) 


8.00 


(0.04) 


0.59 


(0.03) 


14 


17 


-3 


72 


16 


82 


32 


6 


(° 


9) 


9.35 


06538+6355 


6230 


(104) 


8.44 


(0.23) 


0.87 


(0.15) 


14 


67 


-3 


93 


16 


13 


21 


1 


(° 


7) 


9.67 


06556+5920 


18900 


(618) 


7.89 


(0.12) 


0.56 


(0.06) 


10 


70 


-1 


64 


15 


20 


267 


7 


(7 


3) 


7.82 


07029+4406 


11000 


( 46) 


8.02 


(0.05) 


0.61 


(0.03) 


11 


89 


-2 


67 


15 


14 


49 


1 


(1 


4) 


8.68 


07270+1434 


5680 


( 55) 


8.00 


(0.04) 


0.59 


(0.03) 


14 


44 


-3 


82 


16 


58 


30 


7 


(° 


8) 


9.43 


07419+1511 


14370 


(163) 


8.01 


(0.05) 


0.62 


(0.03) 


11 


36 


-2 


20 


15 


18 


75 


7 


(2 


2) 


8.36 


07475+1107 


8160 


( 60) 


8.98 


(0.09) 


1.18 


(0.04) 


14 


62 


-3 


86 


15 


84 


25 


7 


(1 


1) 


9.52 


08126+1737 


16380 


( 80) 


8.09 


(0.05) 


0.67 


(0.03) 


11 


25 


-2 


01 


13 


50 


29 


4 


(° 


8) 


8.24 


09027+3120 


10010 


( 28) 


8.19 


(0.05) 


0.72 


(0.03) 


12 


46 


-2 


94 


14 


75 


35 


3 


(1 


1) 


8.90 


09033+2012 


7200 


( 51) 


8.05 


(0.11) 


0.62 


(0.06) 


13 


49 


-3 


44 


17 


01 


49 


5 


(1 


5) 


9.18 


09106+2156 


5650 


( 59) 


8.00 


(0.04) 


0.59 


(0.03) 


14 


47 


-3 


83 


17 


25 


39 





(1 


0) 


9.44 


09286+1841 


7680 


( 39) 


8.37 


(0.07) 


0.83 


(0.05) 


13 


72 


-3 


51 


16 


51 


37 


6 


(1 


2) 


9.39 


10289+1105 


13910 


(290) 


8.05 


(0.05) 


0.64 


(0.03) 


11 


47 


-2 


28 


15 


19 


64 


9 


(1 


9) 


8.42 


11071+1446 


6620 


( 87) 


8.20 


(0.19) 


0.71 


(0.12) 


14 


04 


-3 


67 


16 


87 


33 





(1 


0) 


9.42 


11401+0112W 


9910 


( 51) 


8.07 


(0.06) 


0.64 


(0.04) 


12 


32 


-2 


89 


14 


60 


49 


2 


(1 


4) 


8.84 


11545+2422 


8770 


(67) 


8.52 


(0.09) 


0.93 


(0.06) 


13 


49 


-3 


39 


15 


26 


26 


3 


(° 


9) 


9.35 


11592+4842 


5900 


( 78) 


8.00 


(0.04) 


0.59 


(0.03) 


14 


26 


-3 


76 


17 


86 


44 


7 


(1 


2) 


9.37 


12155+4630 


6680 


( 78) 


8.12 


(0.17) 


0.66 


(0.11) 


13 


88 


-3 


60 


16 


75 


43 





(1 


2) 


9.33 


12273+3150 


6550 


(67) 


7.97 


(0.16) 


0.57 


(0.09) 


13 


76 


-3 


56 


16 


25 


37 


9 


(1 


1) 


9.23 


13103+1404 


8560 


( 37) 


8.04 


(0.07) 


0.62 


(0.04) 


12 


82 


-3 


12 


16 


08 


49 


7 


(1 


4) 


8.98 


13246+0857 


8110 


( 64) 


7.92 


(0.11) 


0.55 


(0.06) 


12 


86 


-3 


15 


16 


73 


56 


1 


(1 


6) 


8.97 


13349+6945 


14410 


(155) 


7.97 


(0.05) 


0.60 


(0.03) 


11 


30 


-2 


17 


15 


61 


70 


3 


(2 


0) 


8.33 


13521+1053 


6080 


(159) 


8.26 


(0.37) 


0.76 


(0.24) 


14 


50 


-3 


86 


16 


43 


29 


2 


(° 


9) 


9.60 


14106+0245 


5650 


(61) 


8.00 


(0.04) 


0.59 


(0.03) 


14 


47 


-3 


83 


16 


42 


28 


1 


(° 


8) 


9.44 


14149+4336 


6580 


(125) 


8.44 


(0.25) 


0.87 


(0.17) 


14 


43 


-3 


83 


17 


61 


37 





(1 


2) 


9.61 


14553+5655 


15130 


(106) 


8.09 


(0.05) 


0.67 


(0.03) 


11 


38 


-2 


15 


14 


93 


48 


3 


(1 


4) 


8.34 


15206+3903 


8740 


( 44) 


8.04 


(0.07) 


0.62 


(0.05) 


12 


73 


-3 


09 


16 


73 


65 


8 


(1 


9) 


8.96 


15342+0218 


8210 


( 56) 


8.48 


(0.09) 


0.90 


(0.06) 


13 


67 


-3 


47 


16 


43 


32 


9 


(1 


1) 


9.40 


15359+2125 


6560 


(109) 


8.61 


(0.22) 


0.99 


(0.14) 


14 


75 


-3 


96 


16 


87 


38 


2 


(1 


4) 


9.67 


16053+5556 


6810 


( 72) 


7.63 


(0.18) 


0.40 


(0.08) 


13 


16 


-3 


31 


17 


28 


57 


5 


(1 


5) 


9.01 


16096+4735 


13340 


(237) 


8.11 


(0.06) 


0.67 


(0.04) 


11 


63 


-2 


39 


16 


14 


92 


6 


(2 


7) 


8.51 


16171+0530 


14820 


(260) 


7.76 


(0.06) 


0.48 


(0.03) 


10 


95 


-2 


00 


14 


72 


63 


6 


(1 


7) 


8.14 


16264+1938 


6290 


( 93) 


8.39 


(0.21) 


0.84 


(0.14) 


14 


54 


-3 


88 


16 


47 


26 





(° 


8) 


9.63 


16325+0851 


5630 


( 79) 


8.00 


(0.04) 


0.58 


(0.03) 


14 


48 


-3 


84 


15 


05 


13 


2 


(° 


4) 


9.45 


16335+5231 


6670 


( 83) 


8.16 


(0.18) 


0.69 


(0.12) 


13 


95 


-3 


63 


17 


65 


46 


2 


(1 


3) 


9.37 


16546+5742 


9380 


( 63) 


7.76 


(0.09) 


0.47 


(0.05) 


12 


07 


-2 


81 


15 


84 


58 


5 


(1 


5) 


8.73 


17052+0423 


8390 


( 79) 


8.30 


(0.12) 


0.79 


(0.08) 


13 


29 


-3 


32 


15 


95 


35 


5 


(1 


1) 


9.19 


17238+0458 


8450 


( 52) 


8.06 


(0.09) 


0.63 


(0.05) 


12 


90 


-3 


16 


16 


70 


57 


2 


(1 


7) 


9.01 


17283+0211 


8100 


( 47) 


8.13 


(0.08) 


0.67 


(0.05) 


13 


15 


-3 


27 


15 


71 


40 





(1 


2) 


9.10 


17417+2401 


7090 


(61) 


8.33 


(0.12) 


0.80 


(0.08) 


13 


96 


-3 


63 


16 


50 


34 


6 


(1 


1) 


9.46 


18073+0357 


10520 


( 29) 


8.06 


(0.05) 


0.63 


(0.03) 


12 


08 


-2 


77 


15 


08 


31 


6 


(° 


9) 


8.76 



Table 5 — Continued 



PM I 


T B ff (K) 

■*■ err v / 


log 9 


M/Mq 


M v 


loe: L 1 L(7\ 


ya 


D(pc) 


log T 


Notes 


18138+2119 


18100 (112) 


7.97 


(0.05) 


0.60 


(0.03) 


10.90 


-1.77 


15.70 


76.8 


(2.2) 


7.99 




19128+5343 


17670 (100) 


8.32 


(0.05) 


0.82 


(0.03) 


11.48 


-2.03 


13.07 


22.1 


(0.7) 


8.31 




1 Q1 A fi_Ll A OQ 


fiO/i n ?i i ^ 
Dy4U yilo) 


8.60 


(0.21) 


0.98 


(0.13) 


1 A £.1 


— o.oO 


1 T1 
lO.il 


19.4 


(0.7) 






19493+0747 


9410 ( 34) 


8.16 


(0.06) 


0.70 


(0.04) 


12.65 


-3.03 


15.17 


36.3 


(1.1) 


8.95 




20062+5902 


13580 (231) 


7.74 


(0.05) 


0.47 


(0.02) 


11.08 


-2.14 


16.08 


69.3 


(1.8) 


8.26 




20235+7001 


7240 (102) 


8.32 


(0.19) 


0.79 


(0.12) 


13.86 


-3.59 


18.37 


88.0 


'3.0) 


9.42 




21077+0740 


6960 ( 56) 


7.94 


(0.13) 


0.56 


(0.07) 


13.48 


-3.43 


15.91 


34.7 


(1.0) 


9.15 




oi 1 1 7 i ni on 
zlll / -|-UlzU 


100 ( U ^lUz J 


8.06 


(0.05) 


0.65 


(0.03) 


11 1 o 

1 1 . ly 


1 Ofi 

— i.ycs 


10. UU 


67.5 


(1.9) 


q on 

o.ZU 




21134+0727 


6470 ( 66) 


8.20 


(0.15) 


0.72 


(0.10) 


14.14 


-3.71 


15.66 


24.1 


(0.8) 


9.46 




21384+1856 


12610 (225) 


7.65 


(0.08) 


0.43 


(0.04) 


11.10 


-2.22 


16.96 


115.8 


(3.2) 


8.31 


3 


21597+2936 


52110 (839) 


7.66 


(0.08) 


0.56 


(0.03) 


8.54 


+0.34 


14.91 


258.0 


(9.5) 


6.29 


3 


22118+5649 


16790 (166) 


8.15 


(0.05) 


0.70 


(0.03) 


11.30 


-2.01 


12.87 


18.1 


(0.6) 


8.25 




ZZZ f D+l / Oo 


o ( uu ^ yu j 


8.36 


(0.18) 


0.82 


(0.12) 


1 A OA 


— O. I O 


10. o4 


28.8 


(1.0) 


y .oo 




22299+3024 


16310 (157) 


7.48 


(0.05) 


0.37 


(0.02) 


10.38 


-1.66 


15.89 


129.8 


(3.8) 


7.86 


1 


22331+2610 


12020 ( 91) 


7.97 


(0.05) 


0.59 


(0.03) 


11.61 


-2.49 


15.31 


52.4 


(1.4) 


8.55 




23027+4312 


8060 ( 37) 


8.08 


(0.08) 


0.64 


(0.05) 


13.10 


-3.25 


16.85 


79.9 


(2.3) 


9.07 




23160+0559 


15270 (155) 


8.00 


(0.05) 


0.61 


(0.03) 


11.24 


-2.08 


16.93 


138.9 


(4.1) 


8.27 


3 


23229+3358 


14030 (558) 


8.19 


(0.06) 


0.73 


(0.04) 


11.67 


-2.35 


15.68 


101.8 


(3.0) 


8.51 


4 


23234+7255 


7620 ( 47) 


7.99 


(0.09) 


0.59 


(0.05) 


13.19 


-3.30 


16.27 


53.2 


(1.5) 


9.08 




23243+2835 


7130 ( 75) 


8.07 


(0.15) 


0.63 


(0.09) 


13.56 


-3.46 


15.96 


33.4 


(1.0) 


9.20 




23283+3319 


14390 (724) 


8.40 


(0.07) 


0.86 


(0.05) 


11.95 


-2.44 


14.86 


58.9 


(1.9) 


8.62 


4 


23390+5316 


6450 ( 88) 


8.20 


(0.20) 


0.71 


(0.13) 


14.15 


-3.71 


16.79 


32.6 


(0.9) 


9.46 





Note. — (1) Bad GALEX photometry; (2) Binary but only one set of magnitudes; (3) Not in our new selection; (4) DA+dM. 
Estimated from USNO photographic magnitudes. 



